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ABSTRACT 

The orbital angular momentum of a close-orbiting giant planet can be sufficiently large that, if 
transferred to the envelope of the host star during the red giant branch (RGB) evolution, it can spin- 
up the star's rotation to unusually large speeds. This spin-up mechanism is one possible explanation 
for the rapid rotators detected among the population of generally slow-rotating red giant stars. These 
rapid rotators thus comprise a unique stellar sample suitable for searching for signatures of planet 
accretion in the form of unusual stellar abundances due to the dissemination of the accreted planet 
in the stellar envelope. In this study, we look for signatures of replenishment in the Li abundances 
and (to a lesser extent) 12 C/ 13 C, which are both normally lowered during RGB evolution. Accurate 
abundances were measured from high signal-to- noise echelle spectra for samples of both slow and rapid 
rotator red giant stars. We find that the rapid rotators are on average enriched in lithium compared 
to the slow rotators, but both groups of stars have identical distributions of 12 C/ 13 C within our 
measurement precision. Both of these abundance results are consistent with the accretion of planets 
of only a few Jupiter masses. We also explore alternative scenarios for understanding the most Li-rich 
stars in our sample — particularly Li regeneration during various stages of stellar evolution. Finally, 
we find that our stellar samples show non-standard abundances even at early RGB stages, suggesting 
that initial protostellar Li abundances and 12 C/ 13 C may be more variable than originally thought. 
Subject headings: stars: abundances - stars: atmospheres - stars: chemically peculiar - stars: late-type 
- stars: rotation 



1. INTRODUCTION 

Composition studies of red giant stars have yielded 
insightful clues into the internal structure and history 
of low mass stars. Some elemental abundances do not 
change much over the entire lifetime of the star, and 
these abundances can encode information about the 
star's birthplace and the chemical evolution history of 
the gas that became the star. Other abundances, par- 
ticularly those of light elements, are altered by nu- 
clear processes throughout the star's life, and measure- 
ments of these abundances and isotopic ratios have pro- 
vided observational constraints on stellar evolution mod- 
els. Light abundance alterations also provide a means 
of studying planet accretion in red giant stars. The 
abundances of some light element in red giants are sig- 
nificantly depleted from the stars' initial stellar abun- 
dances and, by extension, from the abundances of planets 
that formed from the same protostellar material. There- 
fore, accreted planets can replenish the stellar surface 
abundances of depleted elements. In this work, we fo- 
cus on the lithium abundance, A(Lij3, and the ratio 
of 12 C to 13 C. The initial stellar and planetary abun- 
dances of these elements in solar metallicity stars, as 
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meas ured in the solar system, are A(Li)~ 3.3 dex (see, 
e.g., iBoeseaard et al.lll988tlLodders fe; Feglevll 19981 ) and 
12 C/ 13 C~90 (lAnders fe Grevessdll989lh 

Lithium is destroyed at relatively low temperatures in a 
star by proton capture reactions, and its abundance both 
in the stellar interior and at the stellar surface changes 
considerably throughout the lifetime of most stars. Dur- 
ing the main sequence (MS) phase, Li is preserved only 
in the coolest, outermost regions of the star that com- 
prises only a few percent of the total stellar mass. In the 
stellar interior, the star is essentially devoid of Li. For 
stars more massive than ~ 1.4 Mq, the surface -A(Li) 
stays constant throughout the MS lifetime, while A(Li) 
slowly decreases in lower mass stars. Therefore, by the 
end of the MS, the surface A(Li) in the outer convection 
zone can range from near the original stellar abundance 
to depleted by factors of 100, depending on the star's 
mass. Near the beginning of the post-MS evolutionary 
phases, the star undergoes first dredge-up (FDU), dur- 
ing which the convection zone first deepens into the star 
and then recedes. The Li that had been preserved in the 
outer ~ 1% of the stellar mass on the MS gets mixed 
into the Li-poor inter ior, diluting the surface abundance 
by a factor of ~ 60 (jlbenl Fl967). Similar processes af- 
fect the 12 C/ 13 C ratio. The main chain of the CNO 
cycle reduces 12 C/ 13 C in the stellar interior towards the 
equilibrium value of 3.5, and mixing during FDU will 
lower the overall 12 C/ 13 C of a red gi ant's envelope by 
a factor of ~ 3 ()Sweigart et al.l IT989D . These are the 
"standard model" abundance c hanges, which have been 
obs ervationally v e rified by, e.g., lLambert fe Riesl ()1981[ ) 
and lBrown et all (|1989f ). 

Stars less massive than ~2.3M Q will experience ad- 
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ditional surface abundance changes when the outward 
moving H-burning shell reaches the chemical discon- 
tinuity left behind by the convective envelope dur- 
ing FDU, and this evolutionary stage creates in the 
color-magnitude diagram what is known as the lumi- 
nosity bump or RGB bump. In stars more mas- 
sive than ~2.3¥ 01 the H-burning shell does not reach 
the chemical discontinuity before the star evolves off 
of the red giant branch (RGB). The "erasing" of the 
chemical discontinuity in low mass stars allows for 
additional non-convective mixing methods to further 
decrease the surface A(Li) and 12 C/ 13 C by varying 
amounts. IDenissenkov fc VandenBerel (|2003l ) coined the 
term "canonical extra-mixing" to refer to the prevalent 
processing of light elements beyond standard dilution 
predictions for low-mass, post-RGB bum p giant stars 
that have been observed in open clusters (IGilroyl fl 989; 
IGilrov fc Br owpl ll99lt lLucl[l99l ISmilianic et alH 2009). 
globu l ar clusters dShetrond l2003t iPilachowski et al.l 
120031: iRecio-Blanco fc de Layernvl |2007|) . and fie l d red 
giants dSneden et alj 119861: ICharbonnel et al.l 119981 : 
IGratton et al.ll2000t IKumar fc Reddv l 120091) . The study 
by IDenissenkov fc V andcnBcrg ( 2003) described a two- 
para meter turbulent dif f usion model (following the work 
by iWasserburg et al.1 119951 : iBoothrovd fc Sackmannl 
1999) and ex plored rotation-driven mixing (as previously 
empl oyed by iSweigart fc Mengell979l : ICharbonnel et al.1 
1998) as a physical source of the turbulent diffusion. 
They concluded that most of the non-standard abun- 
dances seen in red giants could be reproduced by a rela- 
tively simple two-parameter diffusion model. 

However, within this framework of standard dilution 
and canonical extra-mixing models, inexplicable light- 
element abundances can remain. Of particular interest 
to the present study are the few percent of red giants 
that ar e "Li-rich," with A( L i) > 1 . 5 dex (Ida Silya et"aLl 
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Gonzalez et alj 2009; Carlb erg et ap2010UlKumar et al 



20111: IMonaco et all |20TT[ IRuchti et all [20TTF 

Lebzelter et al. I 120121) . The Li abundances of these 



stars are larger than that expected from stellar evolution 
and sometimes even exceed the assumed undiluted 
abundances of these stars (e. g.. | Balachandran et all 

1 ICarlberg et al.1 l2010bl: 
Monaco et al.1 120111: IRuchti et al.l 1201 ID . For mas- 
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sive, luminous red giant stars (i.e., M t > 1.5M Q and 
log L/Lq > 4), temperatures at the base of the 
convection envelopes are hot enough for the nucle- 
osynthesis of 7 Li through the Cameron-Fowler chain 
(|Cameron fc FowTerl 119711) . If this freshly synthesized 
Li (or the 7 Be that decays to Li) is quickly trans- 
ported to the cooler regions of the convective envelope, 
then the surface A(Li) can substantially increase 
(|Scalo et al.1 119751) . This mechanism is known as hot 
bottom burning (HBB). However, many of the Li-rich 
giants do not fit this picture because their convective 
envelopes are too cool, even at the hottest depths, 
for L i syn thesis. Both ICharbonnel fc Balachandran 
(2000) and IReddv fc Lambert! (12005) found that many 
of these inexplicably Li-rich giants tend to cluster 
near the luminosity bump. It is thought that the 
erasing of the chemical discontinuity at the bump 



may allow non-convective mixing processes to connect 
the cool convective envelope to hotter depths, where 
Li can be regenerated. Like the HBB mechanism, 
the mixing must be rapid t o replenish the surface 
abund ance. For example, ISackmann fc Boothrovdl 
(1999) used a parameterized "conveyor belt" model 
and could reproduce the abundances of Li-rich giants 
with certain mixing geometries; however, they did 
not provide a physical mechanism for their successful 
models. Rotation was though t to be a likely mechanism 
(IDenissenkov fc Herwigll2004D until lPalacios et al.l (|2006D 
found that a self-consistent model of rotational mixing 
could not generate enough circul ation to build up Li 
in the stellar envelope. Recently, iPalmerini fc Majorca! 
(2010) summarized various mixing mechanisms capable 
of providing the needed conveyor belt and noted that 
the two curren t cont e nders are thermoha l ine m ixing 
(jEggleton et al.1 [2001 ICharbonnel fc Zahnl 120071) and 
magn et ic buoyancy feusso et al.1 12007 \_ Nordhaus et al. 
2008 1 ; iGuandalini et al.l 120091 IPalmerini fc Majorca 
20101) . Because thermohaline mixing is a relatively slow 



process, it may only be an alternative model for the 
canonical extra-mixing that destroys Li. The magnetic 
buoyancy models circulate material fast enough to 
repl enish Li in the s tellar enve lope; however, both 
the | Guandalini et al.1 (|2009l ) and IPalmerini fc Majorca! 
(|2O10l) models predict maximum lithium enrichments 
of A(Li)^2.5 dex, which are an order of magnitude 
smaller than the abundances observed in some of the 
most Li-rich red giants. 

Although there is no consensus on the physical mix- 
ing mechanisms capable of transporting freshly regen- 
erated Li to the stellar envelope, the intrinsically nu- 
clear origin of the Li can be observationally distinguished 
from Li supplied by an accreted planet. First, any ex- 
tra mixing should dredge-up material with low 12 C/ 13 C, 
which will further decrease the su rface 12 C/ 13 C. The "en- 
hance d extra- mixing" models of IDenissenkov fc Herwigl 
(|2004D yield 12 C/ 13 C between 8.2 and 17.4 compared to 
the expected values of 23 without the enhanced mixing. 
Planet accretion, by contrast, should raise both A(Li) 
and 12 C/ 13 C. Second, because only a small fraction of 
stars associated with the RGB bump are Li-rich, then 
either only a small fraction of stars evolving through the 
bump regenerate Li, which is then long-lived in the enve- 
lope, or all stars evolving through the bump phase expe- 
rience a short-lived phase of high Li abunda nces — a "Li 
flash." ICharbonne l & Balachandran (20Q0) considered 
the latter scenario more likely because of the paucity of 
Li-rich giant stars between the luminosity bump and the 
RGB tip. Therefore, Li-rich giants at other evolution- 
ary stages (particularly pre-bump stages) require an al- 
ternative explanation, such as planet accretion. Finally, 
accreted planets should contribute angular momentum 
to the stellar envelope in additio n to altering the stel- 
lar abundances of light elements. ICarlberg et al.1 (2009) 
showed that some known exoplanets are expected to ex- 
perience tidal orbital decay, and the orbital angular mo- 
mentum that is transferred to the stellar envelope is suf- 
ficiently large to measurably increase the host stars' ro- 
tational veloc ities. Because red giant stars are generall y 
slow rotators (|de Medeiros et al.lll996l:lGravlll98lLll982D . 
red giants with enhanced rotation and replenished light- 
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element abundances are candidates for stars that have 
accreted a planetary companion. 

In our project, we select samples of both rapid rotators 
and slow rotators (as a control sample) to look for multi- 
ple signatures of planet accretion simultaneously. In ad- 
dition to a Li enhancement, we are looking for evidence of 
enhancements of 12 C/ 13 C and (in a future paper) refrac- 
tory elements in the atmospheres of rapid rotators, which 
should also be indicative of planet accretion. Interpre- 
tation of these signatures will likely be open for debate 
for any individual star; however, trends that differ signif- 
icantly between the two main samples (defined only by 
their relative i>sini) should make a much stronger case 
to either defend or refute the proposition that some rapid 
rotators gain their angular momentum from a former 
planet. Ours is not the first attempt to correlate Li-rich 
giants and rapid rotation nor the first to invoke planet 
accretion as the explanation for th e correlation. Pla net 
accretion was first put forward by lAlexanderl (Il967fl to 
explai n Li enrichment in giant stars, while [Peterson et al.l 
(1983) were the first to consider plan ets as source s of an - 
gular momentum in evolved stars. Drake ct^l] (|2002f ) 
noted that while both slow and rapid rotators could be 
Li-rich, only 2% of the slow rotators wer e Li-rich com- 
pared to nearly 50% of the rapid rotators. iSiess fc Livid 
(1999) modeled the accretion of sub-stellar companions 
by red giant stars and calculated observational signatures 
of planet accretion in those stars, such as rapid rotation 
and 7 Li enrichment. From the actual occurrence of these 
predicte d observationa l signa tures in the red giant pop- 
ulation, ISiess fc Liviol (|1999D estimated that 4%-8% of 
Sun-like stars host giant planets at orbital radii small 
enough for significant interactions to take place between 
the stars and their planets while the stars are on the 
RGB. Others who have considered planet accretion as the 
solution to b oth Li enhancement and rap i d rotation in 
giants i nclude [W allcrstci n fc Snedenl (119821) .iReddv et al. 
2002) , iCarnev et al.l (|2003t ). and iDenissenkov fc Herwig 
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We describe the selection and observation of our red 
giant sample in Section [2] In Section 02 we detail the 
derivation of the atmospheric properties, rotational ve- 
locity, and abundances. We present our combined abun- 
dance and rotation results in Section 21 where we find Li 
enhancement in the rapid rotators that is consistent with 
planet accretion. We find an unmeasurablc difference in 
12 C/ 13 C between the slow and rapid rotators; however, 
this 12 C/ 13 C result is still consistent with planet accre- 
tion given our measurement uncertainties and the mod- 
est 12 C/ 13 C enhancement expected to accompany the 
observed level of Li enrichment. We divide our sample 
into groups of similar evolutionary stages and find that 
the Li enrichment of the rapid rotators persists within 
these subgroups. In Section [5l we discuss alternative 
internal and external mechanisms that could explain en- 
riched Li and how the enhanced angular momentum of 
the rapid rotators might affect these mechanisms. We 
select a small group of Li rich rapid rotators that are 
the best candidates for planet accretion to test whether 
the rotation and Li enrichment is consistent with planet 
engulfment. Our conclusions are presented in Section [6] 



2. OBSERVATIONS AND DATA REDUCTION 



Our observational sample was built from three dif- 
ferent primary sources. Most of the stars, both rapid 
and slow rotators, come from the sample of Space In- 
terferome try Mission As t rome tric Grid candidates de- 
scribed in lCarlberg et al.l ([201 1L hereafter Cll). Twenty- 
eight candidate rapid rotators were identified from the 
~ 1300 stars studied in Cll. However, because of 
the paucity of the total number of rapid rotators as 
well as the faintness of the Grid Giant Star Sur- 
vey (GGSS) and Tycho stars that constitute the As- 
trometric Grid sample, we supplemented our observ- 
ing list with brighter K giants from the catalog of 
bright F, G, and K giant stars bv lde Medeiros k. Mayor! 
(1999). We included northern hemisphere G9-K9 III 
stars from this catalog that have v smi > 10 kms -1 
and that were not suspected to be spectroscopic bina- 
ries. (Only four stars met all of these selection cri- 
teria and only two were ultimately observed.) Four 
more of our observed ra p id rot ators were selected from 
iGlebocki fc Sta wikowski (2000), which is a compilation 
of rotational velocities for all types of stars. We selected 
bright K II I stars, and their primary v sin i measurement s 
come f rom Strass meier et all (fl9 93). Hen ry et al.l (jl995|) . 
iFekell (|1997f ). and/or iStrassmeier et all (|2000D . A com- 
parison sample of slow rotators was also selected from 
these literature sources. In total, we observed 15 rapid 
rotators and 58 slow rotators. 

In anticipation of measuring the lithium abundance 
and 12 C/ 1 3 C, we adde d a sa mple of known Li-rich gi- 
ants from iDrake et al.l (|2002D and stars with previous 
determinations of 12 C/ 13 C from iDearborn et al.l (I1975D. 
iSmith fc Lambertl ([1985ft , and iFekel &; Balachandranl 
(|1993l) as a control sample to which we can compare 
our abundance determinations. We also found published 
A(Li) for two rapid rotators (HD31993 and HD34198) 
and four slow rotators (Arcturus, HD108255, HD115478, 
and HD116010) after their inclusion in our main sample. 
Finally, we added 13 stars from the literature to act as 
special control samples. Ten of these stars are red gi- 
ants with known planetary companions; they are listed 
in Table [1] together with the associated planetary orbital 
parameters. These stars with planets (SWPs) are all slow 
rotators and are generally expected to show abundances 
similar to those of the other slow rotators (since nei- 
ther sample is expected to have accreted planets). The 
other three stars (HIP35253 HIP3 6896, and HIP81437) 
come from Ma ssarotti et al.l ([20081 hereafter M08), who 
found that these stars' enhanced rotation is most plau- 
sibly explained by planet accretion. These stars form 
an independent sample of stars expected to show planet 
accretion signatures. 

To meet the high signal-to-noise ratio (S/N) and spec- 
tral resolution requirements needed for chemical abun- 
dance work, we obtained spectra with S/N> 100 per 
pixel for our K giant sample using the echelle spectro- 
graphs on the Kitt Peak (KPNO) 4-m Mayall telescope 
and the Apache Point Observatory (APO) 3.5-m tele- 
scope. At KPNO, we chose the 31.6 line mm -1 grat- 
ing and used a 1.5" slit, yielding a spectral resolution of 
R = A/AA w 22,000 when using the long red camera. 
Observations at APO used the default 1.6" slit, yielding 
R ~ 31,500. The observations were taken between 2007 
March and 2010 March. Thorium-argon (Th-Ar) cali- 
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Table 1 

Orbital Parameters for Companions of Stars with Known Planets 



Star Name A4 P sin i P a p e Ref. 

(Mjup) (days) (AU) 



HD 104985 


6.3 


198.2 


0.78 


0.03 


1 


HD 122430 


3.71 


344.95 


1.02 


0.68 


2 


HD 13189 


14 


471.6 


1.85 


0.27 


3 


HD 137759 


9 


511.10 


1.3 


0.712 


4, 5 


HD 177830 


1.5 


410 


1.23 


0.10 


5, 6 


HD 219449 A 


2.9 


182 


2.9 




7 


HD 47536 


5, 7 


43, 2500 


1.61, ■■■ 


0.2, ■■• 


8 


HD 59686 


5.25 


202 


0.911 




9 


HD 73108 


7.1 


269 


0.88 


0.43 


10 


Pollux 


2.7 


589.6 


1.73 


0.02 


11, 12 



spectrum from the magnified image. The target stars are 
fiat-fielded by dividing the extracted target star spectra 
by the normalized flat-field spectrum. 

The extracted, wavelength-calibrated echelle orders 
are then combined and continuum corrected to cre- 
ate normalized, continuous, one-dimensional (ID) spec- 
tra. These ID sp ectra were cross-correlated with the 
iHinkle et al.l (|2000() atlas Arcturus spectrum to measure 
the observed stellar RVs, and these velocities were used 
to shift the wavelength solution to the stellar rest frame. 
The APO spectra span a large wavelength range: there is 
near-continuous coverage between 3200 A and 1 /im. The 
KPNO data, on the other hand, span only 52fO-8250A. 

3. SPECTRAL ANALYSIS 
3.1. Stellar Parameters 

A collection of neutral and singly ionized iron lines can 
be used to constrain the basic stellar parameters of ef- 
fective temperature (T e g), surface gravity (logg), iron 
abundance compared to solar ([Fe/H]), and microtur- 
bulence (£) f ollowi ng the prescription described, e.g., in 
iSmith et al.l (|2001l) . By combining equivalent widths of 
these lines with stellar atmosphere models, one can find 
the model for which all the lines give the same iron abun- 
dance. Generally speaking, Fe I lines spanning a wide 
range of excitation potentials (x or EP) will constrain 
Toff, while a large span of Fe I line strengths will con- 
strain £. Transitions of Fe II are much more sensitive to 
log<? than Fe I, and all of the iron lines will constrain 
[Fe/H] in the form of A(Fe). 

The iron list used to derive the stellar parameters was 
compiled with great care to minimize syste matic errors. 
We began by collecting the i ron li nes from ISmith et al.l 
(|2^0^Jf\ilbright et al.l (|2006l ). and lBizvaev et al.l (|2006(h 
All available log gf values of these sources were tabu- 
lated together with log gf values obtained by querying 
the Vienna Atomic Line Database (VALD; iKupka et al.l 
1999)[] VALD usually returns between one and three 
log gf values per iron line. This initial compilation in- 
cluded 111 Fe I and 17 Fe II lines spanning the wave- 
length range of our spectra. Equivalent widths iW\) of 
these lines were measured using the splot task in IRAF 
for both the solar and Arcturus spectra provided in the 
IHinkle et ail (|2000D spectral atlas CD-ROM. These mea- 
surements were compared to the W\ reported in the 
Smith, Fulbright, and Bizyaev papers. Where our mea- 
surements were discrepant, the lines were remeasured to 
ensure accuracy of our measuring procedure (e.g., veri- 
fying the line identification, checking for overlooked line 
blending). Our line list was large enough that we could 
be discriminating and remove lines that were severely 
blended with neighboring lines or that had a more un- 
certain continuum level. 

We tested our line list by computing solar iron abun- 
dances using the abfind driver of the MOOG stellar 
line analysis prograrrPl (iSnedenl Il973f) together with 
the MARCS (|Gustafsson et all l2008[ isolar atmosphere 
model (T eff = 5777 K, \ogg = 4.40 dex, and A(Fe)= 
7.45). Individual lines that suggested A(¥e) more than 
0.2 dex from the average iron abundance of the entire list 



References. — (1) ISato et al]|2003t (2) ISetiawan et al.H2004 | ; f3l 
Hatzes et all 120051: (4) IFrink et al.l 12003 . disco very; (5) IButler ■ et alj 
200 6|, orbit: (6) IVogt et al.N200q, discovery; (7) [Raghayaii et al.H200a 
f8l| Setiawan et al.M2003l: (9) iMitchell et al.ll2003l : (10) Dolling cr et al.l 
120071 (ll) IHatzes fe Cochranlll99a discovery; (12) IHatzes et al.H2o'oa 
orbit 

bration lamp images, which are used for the wavelength 
calibration, were taken at each telescope position either 
immediately preceding or following each target star's ob- 
servation. For stars requiring very long integration times, 
a Th-Ar comparison image was taken both preceding 
and following the target star observation. Between 2010 
March and 2011 August, we obtained additional low S/N 
echelle spectra of the rapid rotators in an effort to iden- 
tify radial velocity (RV) variations that could signify a 
stellar binary companion as the source of the rapid ro- 
tation. We detail our observations in Table |3l which 
lists each star along with its right ascension (a), declina- 
tion (<5), V magnitude, total exposure time (given as the 
product of the number of exposures and individual ex- 
posure time) , observatory where the data were obtained, 
and the date of observation. The table is organized in 
groups. The slow rotators are presented first followed by 
the rapid rotators, the M08 accretion candidates, SWPs, 
and control stars. Each of the control stars was used as 
a comparison for previous determinations of either A (Li) 
or 12 C/ 13 C. 

The echelle spectra were reduced using standard IRAF 
procedures for overscan correction, bias removal, two- 
dimensional (2D) flat-fielding with quartz lamp images 
(KPNO data only; see below for special flat-fielding pro- 
cedures for the APO data), order extraction, scattered 
light removal, and wavelength calibration. The wave- 
length calibrations of the target star spectra come from 
the associated Th-Ar lamp image(s) taken at the same 
telescope position as the target spectrum. Although 
our targets were observed during bright time, they are 
bright enough that the expected sky flux at full moon is 
less than 1% of the star flux even for the faintest stars 
(V = 12.6). Therefore, we did not perform sky subtrac- 
tion. The APO data required non-standard flat-fielding 
procedures because of how closely spaced the echelle or- 
ders are on the detector. Following the recommendations 
of the APO echelle reduction manual^ we digitally mag- 
nify the 2D images by a factor of four in the spatial di- 
rection before extracting the spectra. We then create an 
average flat-field image, magnify the average flat-field in 
the same manner as the target star images, and extract a 



6 http:/ /www. apo.nmsu.edu/arc35m/Instruments/ARCES/ 
images/echelle_data_reduction_guide.pdf 



7 http:/ /vald. astro. univie.ac.at/~vald/php/vald.php 

8 http:/ /verdi. as. utexas.edu/moog. html 
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Table 2 

Observation Log 



Star 


Q 


S 


V 




Observatory 


DATE-OBS 




(hh:mm:ss.ss) 


(° :':") 


(mag) 


(8) 




(mm/dd/yyyy) 


Slow Rotators 


Arcturus 


i/i i c on Tn 


1 1 n . 1 n c t 

+ iy: 10:57 


0.0 


1 x 0.5 


KPNO 


ni /1 /onno 

01/12/ 2008 


Arcturus 


i /I i c on Tn 

14:15:39. 7U 


1 in in c t 

+ iy: 10:57 


0.0 


1 x (5 + 2) 


KPNO 


no /nc /onnT 

03/06/ 2007 


^nonn i nn o o 

00300+00.29 


no no n i e i 

03:02:04.51 


j nn n 1 nn 

+00:01:00 


12.6 


c ^ ✓ 1 nnn 

b X 1900 


KPNO 


ni /1 1 A /onno 

01/13-14/ 2008 


pri'ii n i cc coon 
O03 iy +56. 5830 


no . o o . oo i n 


1 c 6? . n . n /i 
+ob:2y :04 


10.7 


w nnn 

z x yoo 


KPNO 


ni /in /onno 
01/ 10/2008 


O O 1 O 1 e c c o o o 

O031 ( J+56.6888 


no o c o t e c 

03:2b:27.55 


+56:25:55 


10.6 


x s A n 

2 X 840 


KPNO 


ni /1 /onno 

01/12/ 2008 


/~\ r\ a co i nn n n 

OU45o+UU.9U 


n /i ec oo nn 

04:5b:22.09 


1 nn 1 nn 

+00: 12:00 


12.5 


c \ ✓ icon 

b X lb20 


KPNO 


ni /1 1 /onno 

01/11/ 2008 


(~\ nC On t C C CI 

O063y+56.617y 


n/ 1 a o co cc 

06:43:52. b5 


1 c c nn /in 

+56:0y:4y 


10.4 


w nnn 

2 x yoo 


APO 


11 /1 c /onnT 

11/15/ 2007 


nc r o \ i c ceo 

O0653+16.552 


nc cc cn no 

06:55:50.93 


1 1 c on co 

+ lb:3y:52 


12.3 


5 x 1600 


KPNO 


ni /1 A /onno 

01/14/ 2008 


f~\ C\€£C A 1 1 C O O C 

O0b54+16.235 


nc cc ci /ii 

06:5b:51.41 


+ 16:49:15 


12.4 


5 x 1750 


KPNO 


ni /1 /onno 

01/12/ 2008 


(~~y no a r\ \ e C coon 

00840+56. 583y 


no a c no Tn 

08:4b:08.70 


1 C C 1 A OO 

+56: 14:38 


10.6 


^ / ionn 

2 X 1200 


APO 


no /01 /onnT 

03/31/ 2007 


(~\ n O A r\ \ e c ni on 

O0840+56.9122 


no /ic oo 7i 

08:45:38.71 


1 c c nT c 

+56:07:36 


10.5 


^ / ionn 

2 X 1200 


APO 


no /on /onnT 

03/30/ 2007 


f~* nnnn ne n i 

ooyoy-05.211 


nn. i o./c ot 
U l J: 12:4b. 8 ( 


n C . C T. /I T 

— 05:5 ( :4/ 


11.9 


A w iicn 
4 X 14oU 


KPNO 


n 1 /in/ onno 
Ul/ 1U/Z008 


nn i o n e it 

ooy12-05.11 


nn i c nc nT 

Uy:lo:U5.y7 


—05:58:12 


12.2 


A ^ ✓ ionn 

4 X 1800 


KPNO 


ni /1 /onno 

01/12/ 2008 


^~*nnoe n e i co 

OU9oo-Uo. 152 


nn oo nc oo 

Uy:o8:Ub.22 


n c . a n 10 

— 05:49:18 


12.5 


c ^ / 1 c c n 

6 x lboU 


KPNO 


ni /1 A /onno 

01/ 14/ 2008 


i nc o i nn i c 
OlUOo+UU.lu 


in.ee. /It ic 
10:55:4 1 . lb 


nn. nT.no 
—00:0 f :02 


10.9 


w 1 nnn 
z X 1ZUU 


APO 


no /01 /onnT 
03/31/ zOO / 


riling n r ci 
OllZ4-U0.0l 


1 1 .07,00 TO 

ll:z ( :zo. i 3 


—05:58:53 


12.3 


e w 1 Tnn 
OXI f 00 


KPNO 


n 1 / 1 1 /onno 
01/ 11/Z008 


/~1 i i o t ii c n 

01127-11.60 


ii onoi r* \ 

11:29:21.64 


— 11:24:13 


11.0 


x ✓ ionn 

2 X 1200 


KPNO 


ni /10 /onno 

01/ 12/ 2008 


1 i i o n i O Ct Ct A ~\ A 

01130+39.9414 


ii on i/i on 

11:32:14.20 


1 on n c .ir 

+39:05:4b 


9 . 7 


1 x ✓ ■! nnn 

1 X 4000 


KPNO 


no /nc /onnT 

03/05/ 2007 


/~1 i onn i c T o o o o 

O1200+b7.3882 


ii en on in 

11:59:20.10 


1 C T OC O A 

+67:26:24 


10.7 


x ✓ 1 e nn 

2 X 1500 


KPNO 


no /n/ - * /onnT 

03/06/ 2007 


G 1240+56. 8464 


12:45:02.80 


+55:52:25 


10.3 


1 X 1500 


KPNO 


03/06/2007 


(~ i 1001 i nn 10 

01331+00. 13 


io oo oe eT 

1 3: 33:25.57 


nn no /i c 

—00:03:4b 


10.8 


x ✓ 1 nc n 

2 X 1060 


KPNO 


ni /10 /onno 

01/ 13/ 2008 


/~1 1 in i o o /i c o e 

01421+28.4b25 


i a o -1 in in 

14:24: iy . 10 


1 00 n . nT 

+28:03:07 


10.5 


x ✓ ionn 

2 X 1200 


APO 


no /on /onnT 

03/30/ 2007 


i i e e i i on n .1 e c 

01551 +22. 945b 


i e e a o o Tn 

15:54:22.70 


1 O O O A IO 

+22:24: 18 


10. 1 


1 x 1050 


KPNO 


no /no /onnT 

03/08/ 2007 


/~1 i c .1 n i c/? coot 

Olb40+56.6327 


ic /in .io io 

lb:42:4o. lo 


1 c c nn c t 

+5b:09:57 


10.3 


1 ^ ✓ ionn 

1 X 1200 


KPNO 


no /no /onnT 

03/08/ 2007 


/~i ionn i ci i o n "7 

01800+61. 12976 


io ni ie nc 

18:01: 15.05 


1 C 1 A O CO 

+61:43:52 


10.5 


x ✓ ionn 

2 X 1200 


APO 


nc /no /onnT 

05/02/ 2007 


i^" 1 1 onn i ci i oa7c 
O1800+ bl. 129 to 


io.ni.ic nc 
lo:01: 15.05 


1 C 1 . /I . c 

+Dl:4o:52 


10.5 


1 onn 
Z X 1ZUU 


APO 


ne /on /onnT 
05/ JO/200 / 


s~ i i n o c i ci i -1 o c n 

0193b+bl. 143b9 


ino/i ziooe 

iy:o4:4J.25 


1 c 1 . A n _ a 

+61:49:48 


10.3 


2 x 960 


APO 


ne f r\ a /onnT 

05/04/ 2007 


nnonn i o-icc 

O2z00+5b.34ob 


o i . cn. n c nn 
Zl:o9:05.y0 


1 KG. A A . A Q 

+ 0b:44:4o 


10.6 


Owi onn 
Z X 1ZUU 


KPNO 


n 1 /in/ onno 
01/ 10/2008 


un i noon k 
rl IJlUoZZo 


i o . o c . c n nn 
Iz:zo:o0.y0 


1 on . n 1 . nT 


5 . 


1 x 300 


KPNO 


no /nc /onnT 
03/06/200/ 


i r r "\ -| nnT /I o 

HU1U9742 


io oc co on 

12: 36:58. oO 


1 1 T n c 00 

+ 17:05:22 


5.7 


1 x 300 


KPNO 


no /no /onnT 

03/08/ 2007 


TT T~\ 1 1 r 

HD115478 


io it ic cn 

13:17:15. bO 


1 10 An 00 

+ 13:40:33 


5.3 


1 x 200 


KPNO 


no /nT /onnT 

03/07/ 2007 


un i icnm 
ri.DJ.loU J. U 


i o . on. i n nn 
lo:20: iy .00 


1 -1n.nn.n0 
+4U:U9:02 


5 . 6 


1 x 250 


KPNO 


no /no /onnT 
03/08/200/ 


ttt\i i ooon 

HDiioooy 


i o on no on 

13:39:02.30 


+ 18: 15:55 


6.5 


1 w /mn 1 /ice\ 

1 X (540 + 465 J 


KPNO 


no /ne /onnT 

03/05/ 2007 


ttt\i ni ott 

HD19 1277 


on nc oo oo 

20:05:32.88 


1 c 1 c n a 

+61:59:43 


5.4 


1 X 120 


APO 


ne f r\ a /onnT 

05/04/ 2007 


HD206445 


21:42:10.12 


+01:17:06 


5.7 


1 X 30 


APO 


07/28/2007 


unooi oco 

HD221862 


oo oc nn Ac 

23:35:09.46 


1 c T on . 1 

+67:29:31 


7.2 


2 x 300 


KPNO 


ni /10 /onno 

01/12/ 2008 


tt r\ n c i c n 

HD2blb2 


n a nn i n nn 

04:09: 10.00 


1 in c 00 

+ 19:36:33 


5.5 


1 x 300 


KPNO 


no ! c\r- /onnT 

03/06/ 2007 


rp ni n c n o n o T i 

lycO 195-02087-1 


no in oo cn 

08: 10:33. b0 


1 nn n 1 c T 

+00:01:57 


9 . 3 


1 x 900 


KPNO 


no /nc /onnT 

03/06/ 2007 


T..^rnnc mooT i 

l yclJZUo-UlZo * -1 


no.oo.on 01 
08:2o:3y .31 


1 n a . n . a n 
+04:o0:4y 


9.8 


1 w 1 onn 
1 X 1ZUU 


APO 


no /on /onnT 
03/ JO/200 / 


T.,„noTc nnooT i 

lyc02 /b-0032 /-I 


1 1 .co.no nn 
ll:o8:03.00 


+U4:zl: Iz 


9 . 4 


2 x 600 


APO 


no /on /onnT 
03/ JO/200 / 


T. ,„nQ i n nnooi i 

lycuol9-UUZol-l 


1 a . n a . n c Tn 
14:U4:zo. iU 


1 n /i . no . 
+04:03:2o 


9 . 7 


1 x 900 


KPNO 


no /nc /onnT 
03/06/200/ 


Tyc09 13-0 1248-1 


14:36:30.44 


+ 11:22:11 


10.1 


1 x 1600 


APO 


04/08/2007 


rp ,„nni/i nnC7i i 

lyc09 14-0057 1-1 


i/i /f/ 1 ni on 

14:46:21.20 


+ 11:11:05 


9 . 4 


1 x 1620 


KPNO 


no /nc /onnT 

03/06/ 2007 


Tycl469-0 1108-1 


14:23:16.62 


+ 15:31:38 


9.7 


2 x 660 


APO 


05/02/2007 






4-98- ^9- 


10 1 


9 v 1 son 

^ J\ IOUU 


APO 


1 1 /I o /9D07 


Tycl780-00654-l 


02:43:50.95 


+28:52:58 


10.1 


2 x 500 


KPNO 


01/12/2008 


Tycl890-01314-1 


06:20:21.29 


+28:43:49 


10.0 


2 x 900 


APO 


11/15/2007 


Tycl938-00311-1 


07:59:14.21 


+28:45:54 


9.8 


2 x 660 


APO 


03/30/2007 


Tyc2043-00747-l 


16:24:05.00 


+22:32:30 


9.5 


2 x 900 


KPNO 


03/07/2007 


Tyc2521-01716-1 


10:56:32.22 


+34:39:35 


9.7 


1 x 1200 


APO 


05/30/2007 


Tyc2527-01442-1 


12:10:22.71 


+33:27:00 


9.8 


2 x 500 


KPNO 


01/10/2008 


Tyc3005-00827-l 


10:27:53.44 


+40:16:47 


9.6 


1 x 700 


KPNO 


01/12/2008 


Tyc3013-01489-1 


11:28:35.52 


+39:29:34 


9.9 


2 x 800 


APO 


05/02/2007 


Tyc3027-01042-1 


13:59:16.70 


+37:52:36 


9.9 


1 x 960 


KPNO 


03/08/2007 


Tyc3402-00280-l 


06:53:25.20 


+51:17:55 


9.6 


1 x 1800 


KPNO 


03/06/2007 


Tyc3441-00140-1 


10:15:55.88 


+50:46:44 


9.5 


2 x 480 


APO 


11/15/2007 


Tyc3809-01017-1 


09:14:32.77 


+56:12:38 


9.6 


2 x 600 


APO 


11/15/2007 


Tyc5523-00830-l 


11:51:22.28 


-13:50:18 


11.0 


2 x 1300 


KPNO 


01/12/2008 


Tyc5868-00337-l 


03:05:53.00 


-17:48:48 


9.8 


2 x 500 


KPNO 


01/11/2008 


Tyc5881-01156-1 


03:45:17.54 


-16:29:20 


10.0 


2 x 500 


KPNO 


01/10/2008 


Tyc5981-00414-1 


07:47:19.25 


-16:19:50 


9.6 


1 x 800 


KPNO 


01/10/2008 



were removed. When more than one log gf was available 
for a given line, we chose the log gf value that brought 
the iron abundance closest to the mean. Finally, we re- 
quired that there be no trend of A(Fe) with either the 
lines' excitation potentials or reduced equivalent widths 
(\ogW\/X, or RW). In other words, the "EP-slope" and 
"RW-slope" are both near zero (specifically, within 0.005 
dex eV _1 , standard MOOG parlance). Figure [1] illus- 
trates this zero-slope criterion; output iron abundances 
are plotted on the ordinate against both the excitation 



potential and reduced equivalent width. Both the mean 
A(Fe) and the trend of A(Fe) with the line properties 
are shown. We iteratively removed the lines that gave 
the most deviant A(Fe) until the EP-slope and RW- 
slope were flat and the output A(Fe II) matched the 
output A(Fe I), and this is the solution plotted in Fig- 
ure [TJ The mean solar abundances derived from our 
list is A(Fe I) = 7.53 dex with a standard deviation of 
crpci = 0.08 and A(Fe II) = 7.54dex with cr Fon = 0.06. 
We derive a microturbulence of £ = 1.0 kms -1 . Ten Fe I 
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Table 2 

continued. . . 



Star 


a 


<5 


V 


'"'ex X ^exp 


Observatory 


DATE-OBS 




(hh:mm:ss.ss) 


C :':") 


(mag) 


00 




(mm/dd/yyyy) 



Rapid Rotators 



( ;i isi ).14-SQ 4755 


ns-07-33 nn 

WO.U 1 . O O . W W 


+39: 17:47 


11.9 


2 x 1 80D 

_ s*, xouu 


KPNO 


03 /D6 /?nn7 

WO / WW / aUU 1 


nnsn4— i-^Q 4.7^5 

VJ UOU^r 017.1 1 OO 


no. 07-^^ 00 

UO.Wl .dO.UU 


-I-3Q- 1 7-47 


11 9 


A v 1 ^00 

1 y\ 1 WW 


KPNO 


01/11 /9008 

Ul/ 1 1 1 irWWO 


nn&n4+sQ 47^5 


ns-n7-' : t' : t on 


+39'17"47 


11 9 


1 x 700 


APO 


rn /9Q /9oi 

wo/ iUiu 


G08044-39 4755 


08-07-33 00 


+39:17:47 


11.9 


1 x 700 


APO 


12 /23 /2010 


G0827-16.3424 


08-29-55 02 


— 16:48:13 


9.9 


2 x 500 


KPNO 


01 /I 0/2008 

W 1 / IW/iiWWO 


G0827-16.3424 


08:29:55.02 


— 16:48:13 


f).f) 


1 x 130 


APO 


03/29 /2010 


G0827-16.3424 


08-29-55 02 
wo . o . o o . w ^ 


— 16:48:13 


0.0 


1 x 150 


APO 


10 /22 /2010 


O09984-73 9fsoo 

VJUi7^u~ ' O.^WWW 


09:28:22.03 


■ . w - 00 


10.3 


1 x 1800 


KPNO 


03 /06 /9007 

WO/ WW/ ^UU 1 


009984-73 9600 

VJWy^OH^ ■ u. iUUU 


09:28:22.03 


4-73-09-55 
1 . w y . 


10.3 


2 x 700 


KPNO 


01/11 /9008 
wx/ 11/ ^UUO 


(40Q9S4-7S 9600 


nq.9c.99 OS 


4-7^-00.-^ 
1 . wy . 00 


10 3 


1 x 160 


APO 


0^ /9Q /901 

WO / / iUlU 


(400,9X4-7^ 9600 

Vj W ■} AO 1 O.iUUU 


0Q-9X-99 OS 


1 o.uy. 00 


10 3 


1 x 200 


APO 


1 /99 /901 

1U/^^/ iUlU 


noQ464-oo 4s 

Vj W f 4: U U W . 1 


0Q-4S-41 4S 


_ DO-OQ-OQ 
— uw.uy.u3 


12 4 


V 1 680 

>. J S\ 1 w w 


KPNO 


0^ /07 /9007 


G09464-00 48 


09-48-41 48 


— 00-09-09 
w w . w & . w c 


12.4 


1 x 1100 


APO 


03 /99 /901 

WO/ i J/ iUlU 


noQ.464-00 4s 

Vj UsIU W W . 1 


0Q-4S-41 48 
V , .10.11 .10 


— 00-0Q-0Q 

— UW. \JZ? . Oct 


12 4 


1 x 1100 


APO 


1 9 /9^ /901 

li/iO/ iUlU 


m 21 34-33 1 5558 


12:14:23.35 


4-33-1 1 -45 
1^ .11.10 


11.0 


2 x 1200 


KPNO 


01 /19 /9008 

Wl/ li / iUUO 


HI 21 34-33 1 5558 


12:14:23.35 


4-33-1 1 -45 
1^ .11.10 


11.0 


1 x 300 


APO 


03 /99 /901 

WO / ifC/ / iUlU 


G1 21 34-33 1 5558 

VJ 1 1 O 1 Oi). X O f ) ' .IO 


12:14:23.35 


4-33-1 1 -45 
1 .11.10 


11.0 


1 x 154 


APO 


10 /22 /2010 


G 12 13+33. 15558 


12:14:23.35 


4-33:11:45 


11.0 


1 x 350 


APO 


08/09 /2011 


HD112859 


12:59:03.80 


4-47:09:05 


8.1 


2 x 450 


KPNO 


01/14/2008 


HD112859 


1 2-59-03 80 

li. J^.UO.OU 


4-47-09-05 

TK .UC/.WO 


8.1 


1 x 25 


APO 


03 /99 /901 

WO/ i J/ iUlU 


HD112859 


1 9-59-03 80 


4-47-09-05 

TK .WC/.WO 


8.1 


1 X 25 


APO 


12 /23/2010 


HD112859 


12:59:03.80 


+47:09:05 


8.1 


1 x 360 


APO 


08/09 /2011 


HD31993 


05:00:08.20 


+03:17:12 


7.5 


1 x 780 


KPNO 


03/06/2007 


HD31993 


05:00:08.20 


+03:17:12 


7.5 


1x7 


APO 


10 /22/2010 


HD31993 


05-00-08 90 

UJ.v/U.UO. W 


4-03-1 7-12 

1 W O .1 f . 1 £1 


7.5 


1 x 15 


APO 


12 /23 /2010 


HD33363 


05-18-31 10 

UU.lu.Ol.lU 


4-75-56-49 

' O.OU.I^ 


7.6 


2 x 300 


KPNO 


01 /I 4/2008 

Wl / 11/iiWWO 


HD33363 


05-18-31 10 

Uu.lO.Jl. 1U 


4-75-58-49 

f O.OU.I^ 


7.6 


1 X 20 


APO 


10 /22 /2010 


HD33363 


05-18-31 10 

UU.lu.Ol.lU 


4-75-58-49 

' O.OU.I^ 


7.6 


1 X 20 


APO 


12 /23/2010 


HD34198 


05:14:30.70 


— 26:12:31 


7.0 


1 x 600 


KPNO 


03/07 /2007 


140^41 


O'V 1 4-SO 70 

UO.11.OU. 1 u 


— 96- 1 9-S1 

— u . 1^.01 


7.0 


1 X 10 


APO 


1 /99 /901 

1 W / Ad & I iUlO 


HD34198 


05-14-30 70 
uo.11.ou. 1 u 


— 26:12:31 


7.0 


1 x 15 


APO 


12 /23/2010 


Tvr0' : i47 00769 1 
x y luoi i - ww i \j jL— x 


1 ^-07-49 1 

IO.UI . lil . 1U 


TUO .OO.Ul 


0.8 


9 v 1 900 

_ A. l^UU 


KPNO 


OS /06 /9007 

WO 1 WW/ iUU i 


Tyc0347-00762-l 


1 5-07-49 1 

IO.UI .lii.lU 


4-05-53-01 

TUO . OO - w 1 


0.8 


1 x 120 


APO 


03 /99 /901 

WO/ i J/ iUlU 


Tyc0347-00762- 1 


15:07:42.10 


+05:53:01 


0.8 


1 x 105 


APO 


08/09/2011 


TW0647 009^4 1 
lytuuTii u w^ ji x 


09- c iQ-9fi 71 

U-i . O f . iW. 1 1 


4-1 9-4^-90 

li.lO. m\J 


10 


9 y 600 

i— /K WWW 


KPNO 


m /in /9nos 


Tyc0647-00254-l 


09-59-96 71 


+ 12:45:20 


10.0 


1 X 150 


APO 


10 /22/2010 


Tyc0647-00254-l 


02:59:26.71 


+ 12:45:20 


10.0 


1 x 160 


APO 


12/23/2010 


Tyc0647-00254-l 


02:59:26.71 


+ 12:45:20 


10.0 


1 x 160 


APO 


08/25/2011 


Tyc2185-00133-1 


21:07:09.65 


+28:21:23 


9.8 


2 x 750 


APO 


07/28/2007 


Tyc2185-00133-1 


21:07:09.65 


+28:21:23 


9.8 


1 x 120 


APO 


08/09/2011 


Tyc2185-00133-1 


21:07:09.65 


+28:21:23 


9.8 


1 x 125 


APO 


08/25/2011 


Tyc3340-01195-1 


04:15:10.15 


+51:11:40 


10.0 


2 x 600 


APO 


11/15/2007 


Tyc3340-01195-1 


04:15:10.15 


+51:11:40 


10.0 


1 x 120 


APO 


03/29/2010 


Tyc3340-01195-1 


04:15:10.15 


+51:11:40 


10.0 


1 x 120 


APO 


10/22/2010 


Tyc3340-01195-1 


04:15:10.15 


+51:11:40 


10.0 


1 x 130 


APO 


08/25/2011 


Tyc5904-00513-1 


04:56:02.45 


-17:08:23 


9.5 


2 x 1200 


KPNO 


03/07/2007 


Tyc5904-00513-1 


04:56:02.45 


-17:08:23 


9.5 


2 x 400 


KPNO 


01/10/2008 


Tyc5904-00513-1 


04:56:02.45 


-17:08:23 


9.5 


1 X 90 


APO 


10/22/2010 


Tyc5904-00513-1 


04:56:02.45 


-17:08:23 


9.5 


1 x 95 


APO 


12/23/2010 


Tyc6094-01204-1 


12:05:02.66 


-16:23:42 


9.5 


1 x 720 


KPNO 


01/10/2008 


Tyc6094-01204-1 


12:05:02.66 


-16:23:42 


9.5 


1 x 4000 


KPNO 


03/05/2007 


Tyc6094-01204-1 


12:05:02.66 


-16:23:42 


9.5 


1 x 90 


APO 


03/29/2010 


Tyc6094-01204-1 


12:05:02.66 


-16:23:42 


9.5 


1 x 90 


APO 


12/23/2010 



lines were removed in this process. The remaining lines 
chosen for the final list — 74 Fe I and 13 Fe II — are given 
in Table[3]with the lines' wavelengths, species, excitation 
potential, log gf, and references for the log ^/-values. 

Measuring W\ of the iron lines from the stellar spec- 
tra was accomplished both in an automated and man- 
ual way. A first pass at measuring all of the lines was 
made by interactively running the "Automatic Routine 
for line Equivale nt widths in stellar Spectra" (ARES, 
ISousa et al.l 120071 ). which can automatically fit the con- 
tinuum around each line and fit multiple Gaussians to 
blended lines to measure the W\ of the specified line. A 
visual inspection of each fit is needed to ensure a proper 
treatment of the continuum in line-dense regions. Gen- 
erally, 40%-50% of the ARES-mcasured lines are accept- 
able for stars with v sinz < 11 kms^ 1 . This fraction drops 



to 0%-25% for stars with v sin i > 11km s _1 because of 
severe blending. The remaining lines that are unsuitably 
measured by ARES are manually measured using splot 
in IRAF. The high-resolution Arcturus atlas spectrum 
was used as a guide for identifying the continuum and 
locating individual lines in severely blended regions, par- 
ticularly in the rapid rotators. Not all of these lines are 
measurable in every star. Weaker lines and blended lines 
become increasingly difficult to measure in stars with 
larger vsini values and lower metallicities. In addition, 
non-stellar spectral features such as cosmic rays or tel- 
luric absorption lines can also contaminate some of these 
lines to make them unmeasurable. 

Once the W\ are measured, the stellar atmosphere 
parameters are found with an iterative guess and up- 
date scheme. An initial guess in temperature comes 



Rapid Rotation and Light Element Replenishment 



7 



Table 2 

continued. . 



Star 




5 


V 71 


cx X tcxp' 


Observatory 


DATE-OBS 




(hh:mm:ss.ss) 


(° :':") 


(mag) 


(s) 




(mm/dd/yyyy) 


M08 Accretor Candidates 


HIP35253 


07:17:03.41 


+26:41:22 


6.4 


2 x 200 


APO 


12/02/2008 


HIP36896 


07:35:08.80 


+30:57:39 


5.3 


2 x 10 


APO 


12/02/2008 


HIP81437 


16:38:00.47 


+56:00:56 


5.3 


1 x 30 


APO 


03/29/2010 


Stars with Planets 


HD104985 


12:05:15.12 


+ 76:54:20 


5.9 


1 x 480 


KPNO 


03/06/2007 


HD122430 


14:02:22.78 


-27:25:47 


5.5 


2 x 300 


KPNO 


03/06/2007 


HD13189 


02:09:40.17 


+32:18:59 


7.6 


1 x 800 


APO 


07/28/2007 


HD137759 


15:24:55.77 


+58:57:57 


3.3 1 x 


(240 + 700) 


KPNO 


03/05/2007 


HD177830 


19:05:20.77 


+25:55:16 


7.2 


1 x 800 


APO 


05/02/2007 


HD219449 


23:15:53.49 


-09:05:15 


4.2 


1 x 35 


APO 


07/28/2007 


HD47536 


06:37:47.62 


-32:20:23 


5.3 


2 x 60 


KPNO 


01/11/2008 


HD59686 


07:31:48.40 


+ 17:05:09 


5.5 


1 x 600 


KPNO 


03/05/2007 


HD73108 


08:40:12.82 


+64:19:40 


4.6 


1 x 40 


APO 


03/30/2007 


Pollux 


07:45:18.95 


+28:01:34 


1.2 


1 x (2 + 5) 


APO 


03/30/2007 


Control Stars 


Aldcbaran 


04:35:55.24 


+ 16:30:33 


1.0 


1 X (5 + 20) 


KPNO 


03/05/2007 


HD112127 


12:53:55.70 


+26:46:48 


6.9 


1 x 700 


KPNO 


03/07/2007 


HD127665 


14:31:49.80 


+30:22:17 


3.6 


2 x 60 


APO 


05/30/2007 


HD127665 


14:31:49.80 


+30:22:17 


3.6 


1 x 120 


KPNO 


03/06/2007 


HD163588 


17:53:31.70 


+56:52:22 


3.7 


1 x 30 


APO 


03/30/2007 


HD216228 


22:49:40.82 


+66:12:01 


3.5 


1 x 20 


APO 


05/04/2007 


HD233517 


08:22:46.70 


+53:04:49 


9.7 


1 x 3800 


KPNO 


03/05/2007 


HD33798 


05:15:15.50 


+47:10:15 


6.9 


1 x 2400 


KPNO 


03/05/2007 


HD39853 


05:54:43.60 


-11:46:27 


5.6 


1 x 300 


KPNO 


03/07/2007 


a Stars with 


multiple exposures of 


unequal time arc listed 


as 1 X (*i + t 2 ). These g 


cncrally occurred for the bright stars 


to ensure reaching very 



high S/N without saturating the detector. 
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3 

X (eV) 



■ 4(Fe) = 7.53 
RW-slope=0.001 



-6.5 -6.0 -5.5 -5.0 -4.5 

i°g(f„/A) 

Figure 1. Plot of the output A(Fe) (ordinate) as a function of the 
iron lines' excitation potential (top) and reduced equivalent width 
(bottom). In each panel, the mean A(Fe) for all the lines is plotted 
with a solid line, while the trends of A(Fe) with either \ or W\/\ 
are shown with the dotted lines. 



from each stars' photo metric temperature, using the 
Houdashc lt et al.l (|2000f ) empirical color-temperature re- 
lations for giant stars. The initial guess of the other 
parameters are the same for all stars: logg = 2.0, 
£ = 2.0kms~ 1 , and [Fe/H]=0.0. First, abundances 
are found for a range of T c g and £ around the initial 
guess (holding gravity and metallicity constant) using 
the grid of MARCS stellar atmosphere models. These 
models are defined on a grid of effective temperature 
(3700 < T cf f < 6000, AT cff = 100 or 250 K), surface grav- 
ity (0.0 < logg < 3.5, Alogg = 0.5 dex), and metallic- 
ity (-2.5 < [Fe/H] < +0.25, A[Fe/H] = 0.25 or 0.5 dex). 
When needed, interpolated models are created using the 
interpoLmodeles codes provided with the MARCS mod- 



els0 MOOG fits both the EP-slope and RW-slope for 
each model, and Figure [2] shows an example of how these 
two slopes vary as one runs MOOG on the same line list 
with different T e g and £ but constant \ogg and [Fe/H]. 
A surface is fit to the type of data shown in Figure |2] 
to find the temperature and microturbulence for which 
both slopes are zero. Once these values are found, the 
average A(Fe) given by the Fe I lines is compared to the 
metallicity of the input stellar atmosphere model. If they 
disagree, it is used as the next guess of the metallicity. 
The Fe II lines are used to model a new guess in gravity. 
If the average A(Fe II) is less than the average A(Fe I), 
then a higher surface gravity is guessed in the next itera- 
tion and vice versa for A(Fe II) > A(Fe I). We found that 
updating the log g estimate by the same magnitude as the 
iron abundance discrepancy was a successful scheme. In 
other words, (logcj) ncw « (log.g) old + A(Fe l)-A(Fe II). 
We update T e g and £ only when the EP-slope or RW- 
slope is non-zero by more than 0.005, and we update 
surface gravity only if A(Fe I) and A(Fe II) differ by more 
than 0.03 dex. Our results arc presented in the first 10 
columns of TableHJ which lists the star name, T e s, logg, 
the mean abundance of iron derived from the Fe I lines 
(j4(Fe I)), the standard deviation in that mean (opei), 
the number of lines contributing to the mean (Apd), 
A(Fe II),crFeii, ATpeii, and £. Note that a few stars ap- 
pear in the table more than once because some of the 
stellar parameters are derived independently from more 
than one spectrum of that star. The year of the ob- 
servation is appended to the star name in these cases. 
Stellar parameters that are adopted from one spectrum 
but applied to another (e.g., measured only from a 2007 
spectrum and listed with both the 2007 and 2008 spec- 
tra) are indicated with footnotes. 

9 http:/ /marcs. astro. uu.se/softwarc.php 
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Table 3 

Iron line list 



A Sp 


ocics 


X 


i°g gf 


Reference 


A 


Species 


X 


i°g gf 


Reference 


(A) 




(eV) 


(dex) 




(A) 




(eV) 


(dex) 




5307.361 


Fo I 


1.608 


-2.987 


5, 7 


6127.904 


Fe I 


4.143 


-1.399 


9 


5322.041 


Fe I 


2.279 


-2.803 


5 


6151.617 


Fe i 


2.176 


-3.299 


5, 7 


5412.784 


Fe I 


4.435 


-1.890 


7 


6165.360 


Fe I 


4.142 


-1.550 


5, 7 


5466.988 


Fe I 


3.573 


-2.233 


3, 5 


6173.334 


Fe I 


2.223 


-2.858 


6 


5491.829 


Fe i 


4.186 


-2.188 


1, 5 


6180.203 


Fe i 


2.728 


-2.591 


1 


5522.447 


Fe I 


4.209 


-1.557 


5, 7 


6200.312 


Fe i 


2.609 


-2.437 


5 


5536.583 


Fc i 


2.832 


-3.810 


5 


6232.640 


Fe I 


3.654 


-1.223 


1, 6 


5539.284 


Fc i 


3.642 


-2.660 


5, 7, 9 


6271.282 


Fe i 


3.332 


-2.950 


7 


5560.207 


Fo i 


4.435 


-1.190 


7 


6311.500 


Fe i 


2.832 


-3.230 


7 


5577.031 


Fe I 


5.033 


-1.550 


5, 7 


6322.685 


Fe i 


2.588 


-2.426 


2, 5, 7 


5607.664 


Fo i 


4.154 


-2.270 


5 


6393.601 


Fc i 


2.433 


-1.620 


7 


5611.361 


Fe i 


3.635 


-2.990 


5, 7 


6481.870 


Fe I 


2.279 


-3.113 


8 


5618.633 


Fo i 


4.209 


-1.276 


6 


6518.366 


Fe I 


2.832 


-2.750 


7 


5633.947 


Fe I 


4.991 


-0.270 


5, 7 


6593.871 


Fe I 


2.437 


-2.422 


5, 7 


5635.823 


Fc i 


4.256 


-1.890 


5,7 


6597.561 


Fe I 


4.795 


-0.920 


10 


5636.696 


Fe I 


3.640 


-2.610 


5, 7 


6609.110 


Fe i 


2.559 


-2.692 


5, 7 


5638.262 


Fe i 


4.220 


-0.870 


5, 7 


6699.153 


Fc i 


4.593 


-2.190 


7 


5661.363 


Fo i 


4.284 


-1.736 


5 


6733.151 


Fe i 


4.637 


-1.349 


6 


5691.497 


Fo i 


4.301 


-1.520 


5, 7 


6750.152 


Fe I 


2.424 


-2.621 


2, 5, 7 


5698.023 


Fe i 


3.640 


-2.680 


5, 7 


6820.372 


Fe i 


4.638 


-1.170 


10 


5705.476 


Fe I 


4.301 


-1.600 


7 


6837.020 


Fe i 


4.593 


-1.810 


7 


5712.134 


Fc i 


3.417 


-2.060 


7 


6855.161 


Fe I 


4.559 


-0.741 


3 


5753.122 


Fc i 


4.260 


-0.760 


7 


6858.150 


Fe I 


4.607 


-1.046 


8 


5760.345 


Fe I 


3.642 


-2.490 


5, 7, 9 


6971.937 


Fe I 


3.018 


-3.490 


7 


5778.453 


Fc i 


2.588 


-3.430 


1, 6 


7112.172 


Fc i 


2.990 


-3.090 


7 


5784.660 


Fc i 


3.397 


-2.670 


7 


7189.155 


Fe I 


3.071 


-2.771 


1, 5 


5793.913 


Fe i 


4.220 


-1.829 


8 


7401.658 


Fe I 


4.186 


-1.690 


7 


5807.782 


Fe I 


3.292 


-3.410 


5, 7 


7583.790 


Fei 


3.018 


-1.990 


7 


5809.217 


Fe I 


3.884 


-1.690 


9 


7723.208 


Fo i 


2.279 


-3.617 


2, 5, 7 


5811.917 


Fe i 


4.143 


-2.430 


5 


7941.085 


Fc i 


3.274 


-2.580 


7 


5814.805 


Fe I 


4.283 


-1.970 


5, 7 


5018.440 


Fe II 


2.891 


-1.345 


4 


5837.700 


Fe I 


4.294 


-2.340 


5, 7 


5234.625 


Fe ii 


3.221 


-2.230 


5 


5838.370 


Fc i 


3.943 


— 2.337 


9 


5284.098 


Fc ii 


2.891 


—3.195 


4 


5849 689 


Fe i 


3.695 


-2.990 


5 7 


5^95 559 


Fe II 


3.221 


-3.324 


4 


5853.149 


Fe I 


1.485 


-5.280 


5 


5414.046 


Fc ii 


3.221 


-3.645 


4 


5855.076 


Fe i 


4.608 


-1.478 


1, 5 


5425.247 


Fe ii 


3.199 


-3.390 


4 


5856.084 


Fe i 


4.294 


-1.640 


7, 9 


5991.368 


Fe ii 


3.153 


-3.560 


10 


5861.107 


Fc i 


4.283 


-2.450 


5, 7 


6084.099 


Fe ii 


3.199 


-3.881 


4 


5916.246 


Fe i 


2.453 


-2.832 


6 


6149.246 


Fe II 


3.889 


-2.841 


4 


6024.058 


Fe i 


4.548 


-0.120 


5, 7 


6247.577 


Fo ii 


3.892 


-2.310 


5 


6027.051 


Fe I 


4.076 


-1.089 


5 


6416.921 


Fe ii 


3.891 


-2.680 


10 


6056.005 


Fc i 


4.733 


-0.460 


5 


6432.682 


Fe II 


2.891 


-3.687 


4 


6079.009 


Fc i 


4.652 


-1.120 


5, 7 


6456.383 


Fc ii 


3.903 


-2.185 


4 


6120.250 


Fe I 


0.915 


-5.950 


5, 7 
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The highest surface gravity offered by the MARCS stel- 
lar atmosphere models is 3.5 dex, and one of our stars 
have surface gravities exceeding this upper limit. For 
this star s we switched to using Kurucz plane parallel 
models (|Castelli fc KurucsJ |2001F1 ; a note in Table H 
identifies this star. Because we could not use spherical 
models for all of our stars, we are naturally concerned 
with systematic errors between the stellar parameters de- 
rived with the MARCS spherical models and the Kurucz 
plane-parallel models. For 27 stars in our sample, we had 
derived stellar parameter solutions using plane-parallel 
models previ ous to the spherical m odel derivations pre- 
sented here ([Carlberg et al.l [2010aD . We compare these 
two stellar parameters derivations in Figure |3] The one- 
to-one line is shown for each stellar parameter, and we 
generally find good agreement between the solutions from 
the two different types of models. The only systematic 
difference appears for £, which is slightly lower in the 



10 ODFNEW models 

http: / /kurucz. harvard.edu/grids. html 



available 



at 



spherical models by ~ 0.04 km s -1 , on average. 

To test how well the iron list works on red giant stars, 
we derived stellar para meters for A r cturu s using W\ 
measured both from the iHinkle et al.l (|2000f ) atlas spec- 
trum and our own data. These solutions are tabulated 
along with literature stellar parameters for Arcturus in 
Table [5l This table shows that our line list does well 
constraining T c g, A(Fc), and £. We derive the same 
temperature from the W\ measured in the atlas spec- 
trum and our observed spectrum. Our [Fc/H] differ by 
only 0.03 dex, and £ differ by only 0.04 km s -1 . Compar- 
ing to the literature, we find that our derived T or r devi- 
ate from literature values by < 20 K, our [Fe/H] fall be- 
tween the literature values, and our £ are larger by 0.03- 
0.08 km s -1 . The logg are less well constrained. Our 
two derived values differ by 0.2 dex. Literature values 
range from 1.55 to 1.90 dex, while we derived 1.50 and 
1.70 dex. This larger uncertainty in logg likely stems 
from the fact that log g is primarily constrained by the 
stellar Fe II lines, and there are far fewer Fe II lines in 
our list than Fe I lines. 
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Table 4 

Derived stellar parameters 



Star 


T ff 


log g 


A(Fe i) 


(7f c I 


-NfcI 


A(Fc n) 


(Tfe II 






v sin i 


Error 


N, 


c 




(K) ( 


"111 8 — 2 ) 


(dex) 


(dex) 




(dex) 


(dex) 


(kms -1 ) 


(kms -1 ) 


(kms -1 ) 




(kms -1 ) 


Arcturus_2007 


4300 a 


1.70 a 


6.98 a 


0.11 a 


74 a 


6.98 a 


0.22 a 


12" 


1.64 a 


1.69 


0.61 


5 


5.00 


Arcturus_2008 


4300 


1.70 


6.98 


0.11 


74 


6.98 


0.22 


12 


1.64 


2.36 


0.98 


4 


5.00 


G0300+00.29 


4590 


2.75 


7.52 


0.17 


73 


7.52 


0.23 


12 


1.42 


2.58 


0.84 


5 


4.16 


G0319+56.5830 


5170 


3.10 


7.52 


0.15 


74 


7.51 


0.10 


12 


1.41 


6.92 


0.19 


6 


6.13 


G0319+56.6888 


4790 


2.90 


7.50 


0.15 


74 


7.50 


0.16 


12 


1.55 


2.78 


0.43 


4 


3.12 


G0453+00.90 


4710 


2.80 


7.20 


0.12 


74 


7.21 


0.11 


12 


1.25 


0.66 


0.33 


3 


5.23 


G0639+56.6179 


4660 


2.60 


7.23 


0.10 


74 


7.21 


0.17 


13 


1.15 


1.75 


0.28 


6 


2.80 


G0653+16.552 


4170 


1.40 


7.06 


0.14 


74 


7.05 


0.19 


12 


1.62 


3.25 


0.15 


5 


6.35 


G0827-16.3424 


5090 


3.60 b 


7.36 


0.21 


54 


7.34 


0.15 


7 


2.76 


23.85 


0.42 


5 


3.84 


G0928+73. 2600.2007 


4900 


2.70 


7.26 


0.13 


74 


7.25 


0.16 


12 


1.51 


8.36 


0.30 


5 


5.60 


G0928+73. 2600.2008 


4870 


2.60 


7.28 


0.11 


73 


7.30 


0.11 


12 


1.42 


8.36 a 


0.30 a 


5 a 


5.60 a 


HD31993 


4360° 


2.4 C 


7.62 


0.31 b 


41 


7.43 


0.40 


8 


3.0 C 


30.4 


0.14 


4 


4.59 


Tyc5904-00513-1.2007 


4640 a 


2.30 a 


6.87 a 


0.17 a 


65 a 


6.87" 


0.20 a 


12" 


1.93 a 


14.46 


0.32 


3 


5.09 


Tyc5904-00513-1.2008 


4640 


2.30 


6.87 


0.17 


65 


6.87 


0.20 


12 


1.93 


13.95 


0.32 


3 


7.12 


Tyc6094-01204-1.2007 


4320 


1.90 


7.06 


0.18 


72 


7.06 


0.25 


12 


1.53 


13.71 


0.46 


4 


4.47 


Tyc6094-01204-1.2008 


4320 a 


1.90 a 


7.06 a 


0.18 a 


72 a 


7.06" 


0.25 a 


12" 


1.53 a 


12.78 


0.39 


3 


6.61 


Note. — This table 


is available 


in its entirety in a 


machine-readable form in 


the electronic edition 


of ApJ. 


A portion 


is shown here for 


guidance 



regarding its form and content. 

a Stellar parameter adopted from a previous measurement and not rcmcasurcd from this spectrum 



ments, we find that a Gaussian distribution of a = 0.1 
(i.e., ~ 10% errors) is a reasonable model of our expected 
W\ errors. Therefore, in each iteration of the simulation, 
we vary the equivalent widths of the test star by drawing 
random errors from a Gaussian distribution with a width 
of 10% and rerun our stellar parameter finding software. 
Because this code is relatively time consuming, we only 
ran a small number of iterations. A 10-itcration simula- 
tion for Arcturus revealed mean variations in the derived 
stellar parameters of 80 K in T e g , 0.3dexin log<?, 0.03dcx 
in [Fe/H], and 0.13 km s" 1 in 

3.2. Line Broadening: usini and Q 

All of the stars in our sample have v sin i measure- 
ments in the literature, but these measurements come 
from a variety of sources and from spectra of variable 
quality. Therefore, we re-derive vsini in our entire sam- 
ple. Our high-quality spectra allows us to model the 
shapes of individual lines to measure the broadenings. 
Six relatively isolated lines were chosen from the iron 
lines listed in Tabled A = 5778.45, 6027.05, 6151.62, 
6173.33, 6733.15, and 6750.15 A. For each line, we use 
MOOG's synth driver to generate a synthetic line profile. 
We use the iron abundance measured from that particu- 
lar line, as opposed to the mean iron abundances, to more 
accurately fit the line depth. The broadenings included 
in the profile are the instrumental broadening, macrotur- 
bulence (£), and vsini. The instrumental widths were 
measured from each star's corresponding Th-Ar calibra- 
tion lamp. The widths of five Th-Ar lines in the same 
echelle order as each iron line were averaged together 
to get the instrumental broadening^ The macroturbu- 
lence for each star was estimated from the relations of 

11 In our 2008 KPNO data, the measured instrumental widths 
seemed too large because the broadenings of the slow rotators could 
be reproduced with almost no physical, atmospheric mechanisms. 
One possible reason for this issue is that the seeing conditions were 
better than 1.5" so that the stellar lines are better resolved than 
the ThAr lines. However, that explanation does not account for 
why the 2008 instrumental widths are ~ 40% larger than the 2007 
widths when the same slit was used for both. Within each run, the 



b High log g necessitated the use of Kurucz plane-parallel models. 
c Literature stellar parameters f rom [Castilho et all ]20001) . 
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Figure 2. Illustration of how EP-slope and RW-slope vary with 
temperature and microturbulence for a given log g and [Fe/H] . Solid 
lines connect models of constant T e ff, and dashed lines connect 
models of constant £. The plus sign marks the best-fit T e g and £, 
given in the bottom left corner, as where both slopes are zero. 



Table 5 

Stellar parameters of Arcturus. 



T ff 


logs 


[Fe/H] 


e 


A(Fe) 


Reference 


(K) 


(cm s~ 2 ) 


(dex) 


(kms -1 ) 


(dex) 




4300 


1.70 


-0.55 


1.64 


7.53 


1, 2 


4300 


1.50 


-0.58 


1.68 


7.53 


1, 3 


4300 


1.7 


-0.72 


1.6 


7.50 


4 


4283 


1.55 


-0.50 


1.61 


7.45 


5 


4290 


1.90 


-0.68 






6 



Referenc es. — (1) Derived by us; (2) sp ectrum from this w ork; 
(3) IHinkle et al.h 2000 atlas spectrum: (4) ISmith et al.ll2000l: (5) 
I FuIbright et al.1 12009; (6) IGriffin fc L vnas-Grav 1993, as used in 
iBizvaev et al.112009 . 



We can also constrain our stellar parameter errors us- 
ing Monte Carlo simulations. By comparing the W\ mea- 
sured in our Arcturus spectrum to literature measure- 
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Figure 3. Comparison of the stellar parameters derived using Kurucz atmosphere models and MARCS atmosphere models. The solid 
line is a unity-slope line. 

IHekker fe~ Mclcndezl (|2007f ) based on the stars' temper- 
atures and gravities. Keeping £ and the instrumental 
broadening fixed, we generated spectra with a range of 
v sin i and used a ^-minimization to find the best v sin i. 
Small adjustments to the overall continuum level and ve- 
locity shifts to the wavelength solution were allowed to 
improve the fits. As long as reasonable v smi values were 
measured with the initial £ estimate (i.e., the best fits 
were not v sini=0), then no adjustments to £ were made. 
Otherwise, we reduced £ by using the T a g—^ relationship 
for a dimmer luminosity class, e.g., changing from class 
III to class IV. We do not have the resolution or S/N to 
derive v sin i and £ independently. The average derived 
v sin i value of each star is given in the eleventh column 
of Table @] followed by the standard error in the mean, 
the number of Fe lines comprising the average, and the 
£ of the synthetic spectra. 

To address the question of degeneracy in the v sin i and 
£ derivations, we run the following test. We compute ar- 
tificial spectra at the location of the six iron lines for 
an artificial star with stellar parameters representative 
of our stellar sample, mimicking both the spectral reso- 
lution and sampling of the KPNO data. (The APO data 
are of higher resolution and finer sampling.) Using me- 
dian stellar parameters from our sample, this represen- 
tative star has T eB = 4670 K, log g = 2.55 dex, [Fe/H] = 
-0.21, £ = 1.46 kms -1 , and C = 4.74kms -1 . We are in- 
terested in trying to recover the rotational broadenings 
for four different cases: usini < £, usini w £, usini > £, 

instrumental widths are relatively constant and show only some 
low time frequency variations that are likely due to focus drift. 
These variations are coherent in wavelength. Fortunately, we had 
five stars overlapping our 2007 and 2008 KPNO observing runs. 
By requiring the same v sin i and £ derived from the 2007 data, we 
found that the 2008 KPNO instrumental widths had to be scaled 
by a factor of 0.83. We found no need to change the scaling factor 
with wavelength. 



6749.5 6750.0 6750.5 

A (A) 

Figure 4. Fe lin e at 6750.15 A for the representative star de- 
scribed in Section 13.21 for v sin i of 2, 5, 8, and 20 kms -1 (solid, 
dashed, dot-dashed, and triple-dot-dashed lines, respectively). For 
reference, the approximate instrumental profile (dotted line, a 
Gaussian with FWHM=0.2 A) is superimposed at the central wave- 
length of the Fe line. 

and usini S> C- These cases are modeled with wsini of 
2, 5, 8, and 20 kms -1 , respectively, and in Figure 0] we 
plot the profiles of the 6750.15 A iron line. For refer- 
ence, we have also plotted a Gaussian with a full width 
at half-maximum (FWHM) of 0.2 A. 

For each vsini, we test how the assumed macroturbu- 
lence, Cguess, affects the measured iisini using the same 
v sin i fitting program for the following values of C gU e S s : 
4.74, 2.82, 5.15, and 7.16kms- 1 . The first value is the 
true value of the test star, and we expect to be able to 
recover the correct v sku in these cases. The lat ter three 
values come from the IHekker fe M clcndcz (|2007f ) T-£ re- 
lationships for a T e jf= 4670 K star in luminosity class 
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Figure 5. Difference between the measured vs'mi and model 
v sin i as a function of the model v sin i of a "representative" star. 
Different symbols denote different values of the assumed macrotur- 
bulence, Cguess: the true C (circles), and the expected £ from the 
Hekkcr & Mclcndcz (2007) T — £ relationship for luminosity class 
IV (triangles), III (squares), and II (crosses) stars. 

IV, III, and II, respectively. In Figure [3 we illustrate 
the results of the experiment by plotting the difference 
between the measured and modeled v sin i as a function 
of the model vsini. The different points indicate the dif- 
ferent C's used in the fits. Note that there is no result 
for the vsini = 2kms~ 1 model for Cgucss = 7.164 km s . 
Such a high £ created synthetic lines too broad to fit the 
test lines even in the absence of any rotation. The bar 
in the upper left shows the typical measurement error in 
the fit vsini, which is 0.5kms _1 . As expected, when the 
assumed £ matches the true value of 4. 74 km s -1 , all of 
the model v sin i's are recovered within the error bars. As 
is also expected, the discrepancy between the true vsini 
and the fitted vsini is more dramatic for slow rotators. 
Both rapid rotators in our sample had their v sin i's re- 
covered to within 1.5 km s" 1 across the full range of pos- 
sible Cf 



guess ' 



The slow rotators, on the other hand, may 
have large fractional errors in their derived v sin i val- 
ues if the assumed C is incorrect. However, these larger 
uncertainties merely affirm their classification as "slow 
rotators" and have no other bearing on our analysis. 

In Figure |S1 we compare the v sin i derived here to 
those derived in the literature sources that we used to 
build our sample. All of the comparison vsini val- 
ues for the GGSS and Tycho stars (the fainter sample) 
come from the low S/N spectra described in Cll and 
are shown in the top panel. The brighter stars, plot- 
ted i n the bottom panel , have c omparison measurement s 
fromlde Medeiros et al.l (12000ft . iMassarotti et ail (T2 008) , 
Drake et all (120021 ) ■ iGlebocki fc Stawikowskil (|2000ft , and 
Uesugi fc Fukudal ()1982f ) . In each panel, two comparison 
lines are overplotted. One line represents Avsini= 0, 
and the other indicates the mean difference between the 
vs'mi in the literature and those derived here. We tend 
to measure smaller v sin i compared to literature in both 
the faint and bright stellar samples. The exce ption is the 
subset of stars from de Medeiros et all ((2000), for which 
we measure larger vsini by ~ 0.8 kms' 1 . The most 
likely explanation for the systematically lower v sin i mea- 
sured with our analysis is an overestimation of either the 
instrumental or (more likely) the macroturbulent broad- 



enings. The larger systematic difference seen in the com- 
parison to the Cll data may be due to blended lines in- 
flating the measured line widths in the cross-correlation 
analysis of Cll. In the low- vsini regime (< 5kms _1 ), 
there is larger scatter in Avsini, and trends appear for 
stars with v sin i at the lower limit of the literature sur- 
veys. The higher v sini regime shows less scatter, though 
still a systematic offset. The literature suggests that the 
two most rapid rotators in our sample have even higher 
vsini than what is measured here. Considering the fact 
that the velocity resolution of our data is >9kms _1 , we 
consider the discrepancy between these v sin i measure- 
ments and our comparison samples to be acceptably low. 
In particular, the vsini discrepancy will only affect the 
classification of the stars as "slow" or "rapid" rotators. 
It should have a negligible effect on our abundance mea- 
surements. Measurements of W\ are relatively indepen- 
dent of the line shapes, and spectral synthesis will yield 
correct abundance results in the presence of vsini, £, 
and instrumental broadening errors as long as the total 
broadening is well-represented. 

3.2.1. Defining "Rapid Rotation" 

We choose Skms" 1 as the cutoff vsini between rapid 
and slow rotators. As discussed previously in Cll, the 
v sin i value selected to separate slow rotators from rapid 
rotators is arbitrary, and values between 5 and lOkrns -1 
have been used previously in the literature. Our choice 
of 8 km s _1 is chosen partly because we are confident that 
we can truly distinguish stars with v sin i> 8 km s -1 from 
slow rotators (which was not true for our measurements 
in Cll) and partly for convenience (choosing 8kms _1 
over 7, 9 or lOkrns -1 minimizes the number of "almost 
rapid rotators"). This convention has already been ap- 
plied to the categorization of our targets in Table |31 

3.2.2. Loss of the Independent Planet Accretor Sample 

In our observational sample, the three stars with 
Hipparcos designations (HIP35253, HIP36896, and 
HIP81437) were selected from the v sin i and RV study of 
M08. These three stars have enhanced rotation (between 
7 and lOkms -1 in the M08 study), and the authors ar- 
gued that the location of these stars on the HR diagram 
made the accretion of planets the most likely explanation 
for the high rotational velocities. With reported vsini 
values of 9.9kms- 1 (HIP35253), S^kms" 1 (HIP36896), 
and 7.7 km s _1 (HIP81437) the former two stars would be 
classified as rapid rotators in this study, while the latter 
star just barely misses our cutoff. Therefore, we antici- 
pated that this subsample would be an important control 
sample because their "suspected planet accretor" status 
was made by an independent group. Unfortunately, we 
do not reproduce the v sin i values reported by M08 for 
these three stars. We measure < lkms" 1 (HIP35253), 
S.Skms" 1 (HIP36896), and4.2kms~ 1 (HIP81437). Our 
assumed £ was larger than that of M08 for the two lat- 
ter stars, so we re-derived vsini using their values. Al- 
though this analysis increases the value of our v sin i mea- 
surements, they are still smaller than the M08 values by 
2-3kms _1 and fall short of our rapid rotator cutoff. Be- 
cause we cannot reconcile our v sin i measurements with 
theirs, we simply add these stars to the slow rotator sam- 
ple, which now numbers 61 stars. 



12 



Carlberg et al. 



<i 



15 
10 
5 


15 
10 
5 




avg. An sin i = 2.52 
Carlberg et al. (201 1) 




t\ X - 

ffl + 



avg. An sin i = 1 .76 
de Medeiros & Mayor (2000) 
Massarotfi et al. (2008) 
Drake et al. (2002) 
Glebocki et al. (2000) 
Uesugi & Fukuda (1982) 



_□ 

■O- 



O 
+ 
A 
□ 

X 







10 15 20 

v sin i (km s~ , this work) 



25 



Figure 6. Comparison of the vsini derived in this work to those derived in the literature. T he top panel show s the fainter targets 
selected from Cll (fille d circles), while the b o ttom panel show s the brighter targets selected f rom[dc Medeiros ct al. (2000, ope n circ les), 
IMassarotti et al.l H2008I . crosses') . [Drake et alj (|2002l . triangles). IGlebocki fc StawikowsEI (2000, squares), and lUesugil?F ukuda (1983, *)■ 
In each panel, the lines indicate A«sini= (dotted) and the actual mean A vsini of the samples (dashed). 



3.2.3. Estimating True Rotational Velocity 

Because v sin i depends on the inclination angle of the 
stellar rotational axis to the line of sight, we are in- 
terested in estimating the fraction of slow rotators that 
might be rapid rotators seen nearly pole-on. To test this, 
we computed random inclinations of the rotation axis to 
an observer's line-of-sight for one billion test stars. We 
can then find the fraction of these stars that are true 
rapid rotators (v TO t > 8kms _1 ) given the random incli- 
nation and fsini. This fraction can also be thought of 
as the probability (-Prr) that a star of a given vsini is a 
rapid rotator. We calculate this probability for a range 
of usini between 1 and 5kms _1 , and find that the re- 
lationship between Prr and the observed v sin i is well 
defined by a second-order polynomial of the form: 

Prr = 0.0070 -0.0087(w sin i) + 0.0073(v sin if. (1) 

The expected fraction of slow rotators that are truly 
rapid rotators in our entire sample is simply the sum 
of the probabilities for each slow rotator over the total 
number of slow rotators, which yields 4.5%. Therefore, 
our control sample likely has two or three unidentified 
rapid rotators. 

3.3. Radial Velocity Variability 

One concern with interpreting enhanced rotational ve- 
locity as evidence of planet accretion is the possibility 
that the star is instead merely rotating synchronously 
with a stellar companion. RV monitoring can determine 
whether such a stellar companion is present. There- 
fore, we collected additional high-resolution spectra of 
the rapid rotators to look for RV variability. In Tabled 
we list the star, date of observation, heliocentric RV, the 
uncertainty in the RV (err_RV) , and two measures of the 



likelihood that the star is a single or binary star based 
on the RVs. The first of these two measures, Prv, is the 
probability (from the yf distribution) that the suite of 
measured RVs are consistent with the model of a single, 
constant RV (equal to the weighted mean of the RVs). 
The second measure, /cb, indicates the likelihood of 
measuring Prv for a close-orbiting binary system, given 
the observing cadence and RV precision of each of the 
target stars. To calculate /cb, we generate a popula- 
tion of 50,000 binary systems, drawing random periods 
and mass ratios (g) fr om the Gaussian probability dis- 
tributions described in lDuauennov fc Mavorl (|199ll ). We 
adopted a primary star mass of 2 Mq, and the orbital 
parameters of eccentricity (e), longitude of periastron, 
and time of periastron passage were drawn from uniform 
probability distributions. Inclination angles were drawn 
from an isotropic distribution. We chose a circularization 
period of 100 days (setting e = for shorter periods). 
Because our simulated primary stars are giants, we dis- 
carded simulated systems that yielded binary separations 
smaller than 5 Rq . We define "close orbiting" systems as 
those with periods less than 500 days. These are systems 
for which a 100 Rq primary star that is co-rotating with 
its companion will have w ro t = 10 kms -1 . Each of the 
simulated binary systems is then "observed" using the 
exact same observing dates as the target star, and ran- 
dom noise is added to the RV measurement. These simu- 
lated RVs are then processed identically to the real data 
to compute Prv for each binary simulation. We define 
/cb as the fraction of close binary stars that a have sim- 
ulated Prv within 0.05 of the measured Prv- As an ex- 
ample, in TablelUthe measured P RV for G0827-16.3424 
is 0.06 (or 6%) and /cb = 7%. The interpretation in 
this example is that 7% of the simulated close binary 
systems will have simulated Prv between 0.01 and 0.11 
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Table 6 

Radial Velocity 



Star 


JD-2450000 


RV 


err_RV 


P Rv a 


f b 




(days) 


(km s 1 ) 


(km s — 1 ) 






G0804+39.4755 


4165.660 


12.7 


2.0 


0% 


92% 




4476.817 


35.0 


2.0 








5284.658 


22.3 


0.9 








5553.946 


47.1 


0.9 






G0827-16.3424 


4475.846 


16.1 


2.0 


6% 


7% 




5284.615 


16.7 


1.8 








5492.033 


21.0 


1.4 






G0928+73.2600 


4165.730 


37.9 


2.0 


41% 


0.8% 




4476.894 


37.8 


2.0 








5284.698 


35.9 


0.9 








5492.020 


34.9 


1.0 






G0946+00.48 


4166.661 


44.9 


2.0 


19% 


2% 




5284.679 


40.8 


1.7 








5553.924 


40.7 


1.5 






G1213+33. 15558 


4478.034 


41.4 


2.0 


0% 


91% 




5284.713 


39.7 


0.8 








5492.025 


20.1 


3.8 








5782.606 


40.1 


1.2 






HD112859 


4480.008 


65.6 


2.0 


0% 


89% 




5284.741 


-21.7 


1.8 








5553.909 


58.6 


1.4 








5782.617 


-14.7 


2.1 






HD31993 


4165.611 


17.0 


2.0 


34% 


2% 




5491.925 


14.2 


2.4 








5553.854 


13.2 


1.7 






HD33363 


4479.707 


-57.4 


2.0 


0% 


89% 




5491.990 


-39.9 


1.3 








5553.886 


-37.9 


0.8 






HD34198 


4166.596 


5.4 


2.0 


89% 


0.7% 




5491.932 


5.0 


1.6 








5491.934 


3.7 


1.4 








5553.817 


4.8 


1.4 






Tyc0347-00762-l 


4165.986 


-25.9 


2.0 


0% 


84% 




5284.792 


-29.4 


1.8 








5782.643 


-13.1 


1.8 






Tyc0647-00254-l 


4475.638 


0.1 


2.0 


34% 


0.5% 




5491.909 


-1.1 


0.7 








5553.833 


-1.1 


0.8 








5798.969 


0.8 


O.f) 






Tyc2185-00133-1 


4309.922 


-17.0 


2.0 


71% 


1% 




5782.692 


-15.2 


O.f) 








5798.903 


-15.7 


1.2 






Tyc3340-01195-1 


4419.938 


9.4 


2.0 


0% 


93% 




5284.600 


1.9 


1.3 








5491.974 


7.1 


0.8 








5798.926 


-22.3 


1.4 






Tyc5904-00513-1 


4166.609 


57.8 


2.0 


48% 


0.6% 




4475.712 


57.3 


2.0 








5491.916 


55.2 


1.5 








5553.824 


54.7 


1.3 






Tyc6094-01204-1 


4164.890 


-12.8 


2.0 


0% 


92% 




4476.008 


-28.2 


2.0 








5284.751 


-11.8 


1.1 








5553.955 


-26.1 


1.2 







a Probability that the radial velocities are time independent (RV stable). 

b Fraction of simulated close binary systems (period < 500 days) that have a comparable likelihood of appearing RV stable — Prv (simulated) 
within 0.05 of Prv (actual) — when "observed" with the same cadence and RV precision as the target star. 



(0.06 ±0.05). 

3.3.1. RV Variable Stars and Known Binary Systems 

We find that seven of the rapid rotators have RV vari- 
ability that is inconsistent with the model of a single RV 
(-Prv = 0% in Table [6]). This variability is suggestive 
of orbital motion induced by a stellar companion. One 
RV variable star (Tyc3340-01195-1) is a long-period bi- 
nary star that was known as such before inclusion in 
this study. As was discussed in Cll, the separation 
of the stellar components is large enough that the en- 
hanced rotation cannot be due to synchronous rotation. 
Thus, the rapidly rotating nature of this star is just as 



unusual as that of the other single red giant stars in 
this study. In Figure we compare our RV measure- 
ments of Tyc3340-01195-1 to the RVs pre dicted from the 
known orbital parameters provided by iPourbaix et al.l 
(|200l . and we find that our measurements match the 
predicted curve quite well. Two other RV variable stars 
(HD1 12859 and HD33363) were included in this study 
with the expectation that they were single giant stars, 
but we can confirm their binary nature in the litera- 
ture. HD33363 ap pears in the "9th Ca talog of Spec- 
troscopic Binaries" (jPourbaix et al.l [2004). while the bi- 
nary orbit of HD 112859 was r ecently characterized for 
the first time by iGriffinl (|2009f ) . Our measured RVs are 
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Figure 7. Measured heliocentric RVs (circles) compared to the 
expected velocities from published orbital parameters (solid lines) . 
The abscissa of the top panel is in days, while the bottom two 
panels show orbital phase. 

compared to the expected RVs calculated from the lit- 
erature orbital parameters in the bottom two panels of 
Figure [71 These latter two binary stars have periods on 
the order of only 20 days, and they are likely to be spin- 
ning synchronously with their stellar companions. The 
remaining four stars with variable RVs (G0804+39.4755, 
G1213+33.15558, Tyc0347-00762-l, and Tyc6094-01204- 
1) have not been identified in the literature as binary 
stars before now. The general good agreement between 
the observed and expected RVs in Figure [Jj convinces us 
that this variability is real. At this time, however, we do 
not have enough data to constrain an orbit. The short- 
est time baseline between RV epochs for these four stars 
is 200 days, so we also cannot distinguish whether these 
periods are short enough to be explained by synchronous 
rotation or long enough that the enhanced rotation is still 
unusual (as in the case of Tyc3340-01195-1). 

3.3.2. RV Stable Stars 

The eight remaining rapid rotators have non-zero Pry 
ranging from 6% (possibly RV stable) to 89% (very likely 
RV stable). As we just described in Section [5.3.11 we 
can use /cb to estimate the likelihood that the star 
does in fact have a close-orbiting stellar companion that 
was missed because of the observing cadence. The 
largest /cb of the apparently RV stable stars is 7%, 
for G0827— 16.3424. In other words, if we observed a 
large sample of stars in close binaries with the same 
cadence and observing precision as we have done for 
G0827— 16.3424, we would only expect to measure Prv = 
1-11% in 7% of those systems. The total likelihood of 
this result is even smaller because the "close binary" des- 
ignation applies to only ~ 16% of the binary systems 
(jDuauennov fc Mavorlll991IL and only 14% of t he K gi- 
ants are in binary systems (jFamaev et al.ll2005D . There- 
fore, the total likelihood of G0827- 16.3424 being in a 
close binary may be as low as (0.07)(0.16)(0.14) = 0.2%. 
Because G0827— 16.3424 has the largest value of /cb, we 
conclude that we have observed our targets with suffi- 
cient cadence to detect close binary systems, and stars 
with Prv > are unlikely to have binary companions 
with periods shorter than 500 days. Additional evidence 
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Figure 8. An example fit of a MOOG-generatcd synthetic spec- 
trum (solid line) to the observed spectrum of HD122430 (circles) 
with the location of some spectral features identified. All of the 
labeled atomic and molecular species are neutral. 

of our completeness in detecting binaries is seen in the 
RV variable group (Prv = 0). For all of these stars, more 
than 84% of the simulated close binaries have simulated 
Pav < 5%. This result can be interpreted as a > 84% 
success rate in identifying close binaries with our observ- 
ing cadence. Furthermore, the fact that we detect the 
RV variability of Tyc3340-01195-1 (a long period binary) 
suggests that some of the stars with Prv = may have 
companions that are too distant to cause rapid rotation. 

The two classes of stars relevant to our analysis are the 
rapid rotators that have RV variability that allows the 
possibility of close binary companions (P < 500 days), 
which we refer to as the "close-binary" or CB group) 
and those that are likely single or have only long period 
companions (P > 500 days, which we refer to as the 
"single/long-period" or SLP group). The former group 
contains G0804+39.4755, G1213+33. 15558, HD112859, 
HD33363, Tyc0347-00762-l, and Tyc6094-01204-1, while 
the latter includes G0827- 16.3424, G0928+73.2600, 
G0946+00.48, HD31993, HD34198,Tyc0647-00254-l, 
Tyc2185-00133-1, Tyc3340-01195-1, and Tyc5904-00513- 
1. 

3.4. Lithium 

To measure Li abundances, we use MOOG to gener- 
ate synthetic stellar spectra in a ~ 5 A window spanning 
the Li lines near A = 6708 A. The blending of the Li 
lines with nearby lines of CN, Fe, V, and Ca (see Fig- 
ure [8]) necessitates the use of spectral synthesis over the 
simple equivalent width abundance measurements used 
for the Fe lines. We compiled a line list to compute the 
synthet ic spectra from a vari e ty of sources includin g the 
VALD (IPiskunov et al.lfl995l ). IMandell et al.l (12551 . and 
iGhezzi et al.l ( 2009D . The rest wavelengths, log gf values, 
and excitation potentials for the elements represented in 
this spectral region are listed in Table [7] 

In Figure [8) we show an example synthetic spectral 
fit to the spectrum of HD122430 with the positions of 
major spectral features marked. When fitting, ^4(Li) was 
varied in intervals of 0.01 dex, the smallest abundance 
change specifiable in MOOG. Additionally, we varied the 
abundances of C and N (keeping the ratio C/N constant) 
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Table 7 

Line list for fitting A(Li) 



A 


Species 


X 




(A) 




(eV) 


(dex) 


6 /Ud.o4o 


12 

CN 


3. 130 


— 1.359 


D ( UO.OO ( 


12 /-t]\T 


2. 190 


— 1 .650 


D f Uo.DO f 


12 CN 


0.870 


— 3.001 


6706.733 


12 CN 


0.870 


— 1.807 


6706.844 


12 ,-iTVT 

CN 


1.960 


-2.775 


6706.863 


12 CN 


2.070 


— 1.882 


6706.880 


Fc ii 


5.956 


—4.504 


6706.980 


Si i 


5.954 


— 2.797 


6707.205 


12 CN 


1.970 


— 1.222 


6707.282 


1 2 S~1T.T 

ii! CN 


2.040 


-1.333 


6707.371 


12 CN 


3.050 


-0.522 


6707.433 


Fc i 


4.608 


-2.288 


6707.457 


1 2 ^ii.t 


0.790 


-3.055 


6707.470 


12 CN 


1.880 


-1.451 


6707.473 


Sm ii 


0.933 


-1.477 


6707.518 


V i 


2.743 


— 1.995 


6707.545 


12 CN 


0.960 


— 1.548 


6707.595 


12 CN 


1.890 


— 1.851 


C7fi7 enc 

Diui .oyo 


or i 


4.208 


— 2. 767 


D 1 U 1 .040 




0.960 


— 2.460 


(=17(17 7A(1 


Ce ii 


500 


3 810 


6707 75? 


Sc i 


4.049 


— 2.672 


fi7D7 7^6 

U / \J 1 . 1 ou 


7 Li i 


000 


— 428 


DiUi.i DO 


7 Li i 


n nnn 
u.uuu 


— u.zuu 


6707 771 


Ca i 


5 796 


— 4 015 


K7ri7 sn7 

O ( U / . OU f 


12 CN 


1 210 


— 1 853 


C7H7 Q-l o 




3.600 


— 2.417 


6 ( / .099 


CN 


3.360 


— 3.110 


6707.907 


Li I 


0.000 


— 1.509 


6707.908 


Li I 


0.000 


— 0.807 


6707.919 


Li i 


0.000 


-0.807 


6707.920 


Li i 


0.000 


— 0.807 


6707.930 


12 CN 


1.980 


-1.651 


6707.964 


Ti i 


1.879 


-6.903 


6707.980 


12 CN 


2.390 


-2.027 


6708.023 


Si i 


6.000 


-2.910 


6708.026 


12 CN 


1.980 


-2.031 


6708.094 


V i 


1.218 


-3.113 


6708.147 


12 CN 


1.870 


-1.434 


6708.275 


Ca i 


2.710 


-3.377 


6708.375 


12 CN 


1.979 


-1.097 


6708.499 


12 CN 


1.868 


-1.423 


6708.577 


Fc i 


5.446 


-2.728 


6708.635 


12 CN 


1.870 


-1.584 



as well as Fe, Ca, Si, and V. Small adjustments were 
also allowed in the RV, overall continuum scaling, and 
line broadenings to get the best fits. The quality of the 
fits was estimated using MOOG's plot of the difference 
between the observed and computed spectra. Although 
the total number of free parameters could be numerous, 
we introduced a new free parameter (such as a change 
in the Si abundance or a broadening adjustment) only if 
the spectra could not be well fit without it. For 32 stars 
in our sample, only one of which is a rapid rotator, we 
can only derive upper limits for the Li abundance. 

There are two major sources of uncertainty in the ^4(Li) 
measurements. The first comes from the somewhat sub- 
jective decision of what constitutes the best fit to the 
data. This uncertainty was estimated by seeing how 
much A(Li) could vary before the model was no longer 
a good fit to the data, and we refer to this uncertainty 
as <L4(Li). The second source of uncertainty comes from 
propagating the errors in the stellar parameters to er- 
rors in A(Li). This propagation could not be computed 
analytically because of the complex relationship between 
the stellar structure (as specified in the model stellar at- 



mospheres) and the resulting line profile, and it would 
be excessively time-consuming to remeasure ^4(Li) for 
each star on a grid of different stellar atmosphere mod- 
els in the same manner that we made the initial A(Li) 
measurements. Instead, we used MOOG's abfind driver 
to find the W\ of the Li lines that correspond to our 
measured ^4(Li) from spectral fitting. We then rerun 
MOOG keeping W\ fixed while changing the stellar at- 
mosphere models within the errors of the stellar param- 
eters to measure the variations of A{\A). The T e g and 
£ errors were estimated to be 100 K and 0.5 km s -1 for 
all stars. The errors in A(Fe) and logg are represented 
by ^(Fci) an d <7yi(FciiJ3 from Table 21 respectively. For 
most stars, only the variation in T c g yielded variations in 
A(Li) greater than 0.01 dex (the output abundance pre- 
cision of MOOG). The total formal error for each star, 
e vi(Li)) is the quadrature sum of the propagated stellar 
parameter error and 5A(Li). 

Finally, these Li abundances, which were derived in 
local thermodynamic equilibrium (LTE), were corrected 
for non-LTE (NLT E) departur e s usin g the grid of cor- 
rections supplied in iLind et ail ()2009l ) . The grid points 
include 4000 < T eff < 8000 K with AT off = 500 K, 
1 < logg < 5 with Alogg = 1, -3.0 < [Fe/H] < 
0.0 with A [Fe/H] = 1.0, -0.3 < A(H) < 4.2 dex with 
AA(Li) = 0.3, and £ = 1, 2, and 5kms _1 . We used sim- 
ple linear interpolation to estimate corrections between 
the Lind grid points to match the stellar properties of 
T g, logg, [Fe/H], and A(L1)lte, and we chose to use 
the grid that most closely matched each star's £ (instead 
of interpolating between £) . A number of stars had stel- 
lar parameters outside of the grid, most commonly for 
having A(Li)LTE < —0.3 dex. However, we noticed that 
when all other stellar parameters were held constant, 
the NLTE corrections tended toward a constant value 
for small A(L1)lte- Therefore, we used the corrections 
supplied for A(L1)lte = —0.3 dex for these low lithium 
stars. Similarly, we used the corrections computed at 
[Fe/H]= 0.0 for the 13 stars with super-solar metallicity 
and the corrections computed for log g = 1 for Tyc0913- 
01248-1, which has logg = 0.9 

However, we found that the ILind et al.l ()2009[ ) NLTE 
correction grids were not completely filled over the stel- 
lar parameter ranges described above, which led to a 
problem for 44 stars using the £ = lkms -1 grid. The 
interpolations would have included "empty" NLTE cor- 
rections in multiple dimensions of the grid. The prob- 
lem was rectified for 43 of these problem stars by us- 
ing the £ = 2kms~ 1 grid. Over the stellar parameter 
range covered by our target stars, the NLTE corrections 
in the £ = lkms -1 and £ = 2kms _1 grids (when cor- 
rections existed for both) differed by only 0.008 dex. For 
the remaining problem star, HD177830, we simply re- 
port the LTE abundance. Our final Li measurements are 
given in the first five columns of Table 13 which lists the 
star name, the measured LTE abundance, 6 A (Li), e^(Li), 
and the NLTE-corrected ^4(Li) for each of our program 
stars. Upper limits to A(Li) are indicated with "<" in 
the <L4(Li) column. 

As a check of our Li measurements, we compare our 

12 Recall from Section 13.11 that we found that changing logg 
of the stellar atmosphere models by A log g effected a change in 
A(Fe II) that was ~ - A logg. 
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Table 8 

LTE and NLTE A(Li) and 12 C/ 13 C of the program stars 



Star 


A(Li) L 


<5A(Li) 


<U(Li) 


A(Li) N 


12 C/ 13 C 




(dex) 


(dex) 


(dex) 


(dex) 




Arcturus_2007 


-1.23 


< 


0.17 


-0.93 


7.5+0.3 


Arcturus_2008 


-0.73 


< 


0.17 


-0.43 


6.9 ±0.5 


G0300+00.29 


-0.19 


0.2 


0.25 


+0.05 


40.3 ± 5 


G0319+56.5830 


+0.01 


< 


0.11 


+0.14 


17.7 ± 3.6 


G0319+56.6888 


-0.61 


< 


0.13 


-0.41 


14.0 ± 3 


G0453+00.90 


+0.36 


0.02 


0.13 


+0.55 


21.5 ± 5 


G0639+56.6179 


-0.99 


< 


0.14 


-0.78 


19.0 ± 1.2 


G0653+16.552 


+0.35 


0.01 


0.16 


+0.70 


11 ± 3 


G0654+16.235 


+0.13 


0.01 


0.16 


+0.43 


11.9 ± 1.5 


G0804+39.4755 


+ 1.24 


0.01 


0.13 


+ 1.41 


> 30 


G0928+73.2600 a 


+3.62 


0.05 


0.20 


3.30 


28 + 8 


HD31993 b 


+0.98 


0.1 


0.19 


+ 1.03 


13.6 ± 2 



Note. — This table is available in its entirety in a machine- re ad able 
form in the electronic edition of ApJ. A portion is shown here for guid- 
ance regarding its form and content. 

a See lCarlberg et all d2010bn for details of the A{Li) and 12 C/ 13 C mea- 
surement . 

b Using literature stellar parameters from Castilho et al. (2000). 



Table 9 

Comparison of yl(Li)LTE measurements. 



This Work Literature 



Name 


A(U) 


v sin i 


A(Li) 


v sin i 


HD34198 


0.08 + 0.2 


17.3 


0.4+0.2 


18.7" 


HD33798 


1.66+0.1 


28.0 


1.5+0.2 


30 b 


HD31993 


0.98°, 1.80+0. 2 


30.4 


1.4+0.2 


31.1* 


HD108255 


< -0.49 


0.1 


< +0.6 d 


1.4 b 


HD115478 


-0.36 ± 0.2 


0.0 


-0.4 d 


1.7 b 


HD116010 


< -1.5 


0.5 


< -0.6 d 


1.4 b 


Arcturns 


< -0.73 


2.0 


< -0.8 d 


< 1.0 b 


HIP81437 


-0.03 ± 0.2 


4.2 


< +0.2 d 


8 e 



" IBohm -Vitcnsc I2004T). 

b Ide Medeiros fc Mavorf II1999D . 

c Our derivation using stellar parameters from Castilho et al. (2000). 
d IBrown et al.H1989|). 
" IMassarotti et al.l J2OO8I) . 



measurements to published values for eight stars, as 
listed in Table [H Because the literature A(Li) were de- 
rived under the assumption of LTE, we use our LTE- 
derived values for comparison. The literature sources 
also each provide v sin i so we include a v sin i compar- 
ison in Table [§] as well. Only four stars have literature 
j4(Li) measurements that are not upper limits. Of these, 
all of the Li measurements agree within the quoted un- 
certainties. 

3.5. Carbon Ratio 

The measurement of 12 C/ 13 C comes from fitting a 
small group of lines in the spectral region between 8001 
and 8005 A, where synthetic spectra are again generated 
with MOOG. The line list used for spectral synthesis is 
presented in TableHUl which lists each line's wavelength, 
species, excitation potential, log gf, and the primary ref- 
erences. The CN wave lengths are taken f rom t he labora- 
tory measurements bv lDavis &: Phillipi ([1963) for 12 CN 
and iWvllerl (fl966l ) for 13 CN. The gf values were calcu- 
la ted based on the ab s olute / values from the analysis 
of iSneden fc Lambertl (|1982l ) and the adopted dissocia- 
tion energy is 7.65 eV. The atomic line list was complied 
using VALD (jPiskunov et al.lll995l ). The log gf value of 



the Fe I lines at 8002. 567 A and 8003. 227 A are reduced 
by ~ 0.6 dex from the values found in VALD. Without 
this reduction, the models predict significantly stronger 
absorption at these wavelengths than what is present in 
the observed spectra — both in the atlas solar spectrum 
and our red giant spectra. 

The strengths of the molecular features we are fitting 
depend on the abundances of both C and N, and we al- 
low both of these abundances to chan ge. However, we 
fix th e ratio of C to N to 1.5 (following IGustafsso n et al.l 
l2008f ). which is appropriate for moderate CN processing 
in giant stars. This method allows us to add only one 
free parameter to describe C and N abundances that de- 
viate from scaled-solar. Fitting this red region of the 
spectrum is also complicated by the presence of numer- 
ous telluric features, mainly H2O. No telluric calibration 
stars (such as a rapidly rotating [wsini > 200 km s -1 ] O 
or B stars, which have nearly featureless continua) were 
observed with the target stars; however, the high qual- 
ity atlas telluric spectrum from iHinkle et al.l (j2000) pro- 
vides a good model. Although it is common to "divide 
out" the telluric features that are superimposed on stel- 
lar spectra, we instead chose to add a telluric spectrum 
that is scaled, broadened (to maintain the instrument 
response), and RV shifted (to place it in the stellar rest 
frame) to the synthetic spectrum models. To account 
for the exponentially decreasing transmission of the stel- 
lar spectrum through the atmosphere at wavelengths of 
telluric absorption, i.e., applying Beer's Law, the scaled 
telluric spectrum T(A) is given as T(A) = T (\y, where 
To (A) is the broadened but unsealed telluric spectrum 
and / is the scaling factor, a free parameter. 

Other free parameters in the overall fit include small 
adjustments in the stellar broadening (generally within 
2kms _1 each for the vsini and (), RV (usually less than 
2kms _1 ), iron abundance (up to ± 1 o\4(Fei))> an d con- 
tinuum adjustments. The continuum fitting is the most 
subjective and, in some cases, the most difficult part of 
the fitting procedure because of the high density of stel- 
lar lines. Unlike the ^4(Li) fitting, where a global scaling 
factor was sufficient to get good fits to the continuum, 
we allow linear adjustments to be made. Higher order 
fits were considered, but ultimately it was decided that 
higher order adjustments were more likely to introduce 
errors rather than reduce them. For stars with extreme 
rotation (vsmi > 15kms _1 ), there is no continuum in 
the spectral region being fit. In these cases, we relied 
heavily on the models for guidance. 

A best-case example of measuring 12 C/ 13 C in our stel- 
lar spectra is shown in Figure |H] for G1130+37.9414. In 
the top panel, the three synthetic spectra demonstrate 
the best fit 12 C/ 13 C of 25 and 12 C/ 13 C that are both 
higher and lower than the best fit. Major spectral fea- 
tures are identified. The synthetic spectra include a tel- 
luric contribution, which is plotted separately for refer- 
ence. In the bottom panel of Figure HI we show the best 
fit decomposed into the pure stellar spectrum that is gen- 
erated with MOOG and the atlas telluric spectrum. We 
used the strong telluric line indicated in that panel to 
find the scaling parameter, /, for the telluric spectrum. 

In Table ITTT we compare our measured 12 C/ 13 C (from 
the technique described above) to published values for 
five of our control stars and three of our program stars. 



Rapid Rotation and Light Element Replenishment 



17 



Table 10 

Line list for fitting 12 C/ 1 



A 


Species 


X 


log gf 


Reference 


A 


Species 


X 


l°g 9f Reference 


(A) 




(cV) 


(dcx) 






( A ) 




(cV) 


(dex) 




7990.388 


12 CN 


1.38 


-2.0585 


7, 


8 


8002.367 


13 CN 


1.49 


-1.8327 


7, 9 


7990.790 


12 CN 


1.45 


-1.6234 


7, 


8 


8002.412 


12 CN 


0.18 


-1.4962 


7, 8 


7990.957 


13 CN 


1.37 


-1.6861 


7, 


9 


8002.571 


13 CN 


0.21 


-1.7212 


7, 9 


7991.128 


13 CN 


0.05 


—2.1512 


7, 


9 


8002.576 


Fe i 


4.580 


-2.2400 


2 


7991.128 


13 CN 


0.18 


-1.7496 


7, 


9 


8003.185 


Al i 


4.087 


-1.8791 


4 


7991.583 


12 CN 


1.47 


— 1.6216 


7, 


8 


8003.213 


12 CN 


0.12 


-1.9431 


7, 8 


7991.711 


13 CN 


0.04 


-2.2700 


7, 


9 


8003.227 


Fc i 


5.539 


-2.3889 


2 


7992.126 


13 CN 


0.10 


-1.6364 


7, 


9 


8003.311 


13 CN 


1.34 


-2.0883 


7, 9 


7992.297 


12 CN 


0.09 


—2.0114 


7, 


8 


8003.485 


Ti i 


3.724 


-0.2000 


1 


7992.297 


12 CN 


0.18 


— 1.5143 


7, 


8 


8003.553 


12 CN 


0.31 


-1.6440 


7, 8 


7993.144 


13 CN 


0.21 


— 1.7423 


7, 


9 


8003.910 


12 CN 


0.33 


-1.6478 


7, 8 


7993.600 


Ti i 


1.873 


— 2.4970 




1 


8004.036 


12 CN 


0.06 


-2.9245 


7, 8 


7993.869 


13 CN 


1.49 


— 1.8508 


7, 


9 


8004.550 


13 CN 


0.12 


-1.5918 


7, 9 


7994.018 


12 CN 


1.42 


— 1.6326 


7, 


8 


8004 71 5 


13 CN 


0.07 


-2.0814 


7, 9 


7994 31 9 


13 CN 


0.09 


— 1.6615 


7, 


9 


8004 sol 


13 CN 


0.10 


-1.6144 


7, 9 


/ yyi.ooo 


13 CN 


1 34 


^ 6536 


7, 


9 


OUUO . Zlo 


Zrl 


623 


-2.1901 


4 


7QQ4 fiQA 

/ yy^.oy^ 


12 CN 


11 


^ 9666 


7, 


8 


OUUD.UDO 


13 CN 


1 41 


-1.6517 


7, 9 


l iJiJO .KJ X<J 


12 CN 


16 


— 1 5143 


7, 


8 


snnfi 1 9(=; 

OUUD . -LZO 


13 CN 


24 


-1.7122 


7, 9 


/ yy □ . too 


13 CN 


02 


9 RQQfi 
— z .oyyo 


7, 


9 


snnfi 4^q 


Si i 


6 261 


-1.7231 


5 




12 CN 


06 


— 2 9172 


7, 


8 


snnfi 7n^ 


Fc i 


5 067 


-2.1280 


2 


7qqc; (\An 


12 CN 


29 


— 1 6556 


7, 


8 


ovjuu . y.£tj 


12 CN 


1 60 


-1.7878 


7, 8 




12 CN 


31 


— 1 6615 


7, 


8 


8no7 9i i 


13 CN 


1 48 


-1.2652 


7, 9 


7996.435 


Ti i 


3.337 


0.2660 




1 


8007.242 


Co i 


4.146 


0.1159 


6 


7996.706 


13 CN 


0.19 


— 1.7352 


7, 


9 


8007 ^82 


12 CN 


0.11 


-1.9586 


7, 8 


7996.761 


12 CN 


1.58 


— 1.7986 


7, 


8 


8007.882 


13 CN 


0.03 


-2.8962 


7, 9 


7996.807 


13 CN 


0.04 


— 2.2233 


7, 


9 


8007.904 


13 CN 


0.05 


-2.1415 


7, 9 


7996.816 


Fc i 


4.584 


-2.4860 




2 


8008.387 


Si i 


6.079 


-1.8289 


5 


7997.334 


13 CN 


1.39 


-1.6696 


7, 


9 


8008.493 


12 CN 


0.20 


-1.4802 


7, 8 


7997.800 


13 CN 


0.06 


-2.1152 


7, 


9 


8008.652 


13 CN 


1.36 


-1.6421 


7, 9 


7998.216 


13 CN 


0.11 


— 1.6144 


7, 


9 


8008.737 


13 CN 


0.22 


-1.7077 


7, 9 


7998.312 


12 CN 


1.54 


— 1.2480 


7, 


8 


8009.278 


12 CN 


1.44 


-1.6216 


7, 8 


7998.876 


13 CN 


1.38 


-1.6696 


7, 


9 


8009.703 


12 CN 


1.62 


-1.7932 


7, 8 


7998.944 


Fe i 


4.371 


0.1489 




2 


8010.084 


12 CN 


0.19 


-1.4802 


7, 8 


7999.214 


12 CN 


1.40 


-2.0287 


7 


g 


8010.084 


12 CN 


1.40 


— 9 0900 


7 8 


7999.214 


12 CN 


1.60 


-1.8041 


7, 


8 


8010.468 


13 CN 


0.11 


— 1.5918 


7, 9 


7999.408 


13 CN 


0.09 


-1.6383 


7, 


9 


8011.101 


13 CN 


1.51 


-1.8210 


7, 9 


7999.460 


13 CN 


1.46 


-1.2644 


7, 


9 


8011.159 


13 CN 


0.14 


-1.5702 


7, 9 


7999.465 


13 CN 


0.22 


-1.7282 


7, 


9 


8011.484 


13 CN 


1.63 


-1.4828 


7, 9 


7999.846 


12 CN 


0.10 


-1.9830 


7, 


8 


8011.732 


12 CN 


0.33 


-1.6326 


7, 8 


8000.261 


12 CN 


0.19 


-1.4962 


7, 


8 


8011.836 


13 CN 


1.52 


-1.8268 


7, 9 


8000.316 


12 CN 


1.47 


-1.6091 


7, 


8 


8011.899 


13 CN 


0.08 


-2.0506 


7, 9 


8000.757 


Nd ii 


1.091 


-1.2220 




3 


8011.950 


12 CN 


0.13 


-1.9208 


7, 8 


8001.369 


13 CN 


0.03 


-2.8962 


7, 


9 


8011.950 


12 CN 


1.50 


-1.5952 


7, 8 


8001.524 


12 CN 


1.42 


-1.6253 


7, 


8 


8012.546 


12 CN 


0.35 


-1.6364 


7, 8 


8001.652 


12 CN 


1.48 


-1.6091 


7, 


8 


8012.620 


13 CN 


1.35 


-2.0675 


7, 9 


8002.214 


13 CN 


0.05 


-2.1805 


7, 


9 
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One control star, HD39853, had a strong telluric absorp- 
tion feature fall precisely at 8004.7 A in the stellar rest 
frame, rendering 12 C/ 13 C unmeasurable. For the remain- 
ing stars, Table HT1 demonstrates that we do tend to mea- 
sure similar 12 C/ 13 C to those found in the literature. For 
example, with the exception of Aldcbaran, the ranking of 
stars from low to high 12 C/ 13 C would be the same using 
our measured 12 C/ 13 C or using the literature 12 C/ 13 C. 
However, the measurements often disagree by more than 
the quoted uncertainties, indicating that either we or the 
literature sources are underestimating the uncertainty. 

We add the derived 12 C/ 13 C and uncertainty for our 
stellar sample in the final column of Table [SJ The uncer- 
tainties for some of the stars come from averaging sev- 
eral different attempts at fitting the line that all gave 
equally good results. In these situations, the listed un- 
certainty is the standard deviation in the individual mea- 



surements. In the absence of multiple measurements we 
assign an uncertainty of ±3 for 12 C/ 13 C < 20 and ±5 
for 12 C/ 13 C > 20. These choices come from the median 
uncertainties in the stars with multiple measurements. 



4. RESULTS: EVOLUTION OF ABUNDANCES AND 
ROTATION 

In this paper, we are investigating the hypothesis that 
rapid rotators can be spun up by the accretion of planets 
by searching for the chemical signatures of 4(Li) replen- 
ishment and increased 12 C/ 13 C that should accompany 
the enhanced rotation under the planet accretion sce- 
nario. If the enhanced rotation is not caused by planet 
accretion, then we expect to find that either (1) 4(Li) 
and 12 C/ 13 C will not differ between the rapid rotator and 
slow rotators samples, or (2) if rapid rotation drives addi- 
tional mixing conducive to lithium production, any 4(Li) 
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Figure 9. Top: sample fits to measure 12 C/ 13 C from the observed 
spectrum of G1130+37.9414 (circles). The three synthetic spectra 
have 12 C/ 13 C of 40 (dashed), 25 (solid, best fit), and 15 (dotted). 
Major spectral features are labeled with the contributing atomic or 
molecular species. The weak 13 CN lines near A = 8004.7 A provide 
the leverage for measuring 12 C/ 13 C. The lHinkle et al.l ||2(300T) atlas 
telluric spectrum, broadened to match our instrumental resolution, 
is shown for reference (triple-dot-dashed line, offset vertically for 
clarity). Bottom: redward of the 13 CN feature is a telluric line 
in the observed spectrum (circles) that is suitable for scaling the 
atlas telluric spectrum. The lines show the best-fit 12 C/ 13 C= 25 
model (solid) and the contributions from the MOOG-gencratcd 
stellar spectrum (dashed) and the atlas telluric spectrum (triple- 
dot-dashed). 



Table 11 

Comparison of 12 C/ 13 C measurements 



Name 


12 C/ 13 C 




This Work 


Literature 


Aldcbaran a 


19.7 ± 3 


10 b 


Arcturus 


7.2 ± 0.2 


7° 


HD112127 


34.4 ± 2 


22 d 


HD127665 


11.3 ± 0.5 


15 ± 2° 


HD163588 


25.8 ± 2 


20 ± 2° 


HD216228 a 


11.4 ± 0.4 


16 ± 2° 


HD39853 


Telluric contamination 


7 d 


Pollux 


15.9 ± 3 


16 c 



a Used literature values fo r model atmosphere. 
b ISmith fc Lamberti J1985D . 
c IGilrovl IU989TT 

d IFekel &: Balachandrajil 119931) . 
' IDearborn et al.l U975I) . 



enhancement will be accompanied by lower 12 C/ 13 C. We 
explore these scenarios by comparing both the group 
abundance properties of the slow and rapid rotators (Sec- 
tion 14. 1|) and dividing the samples into common stellar 
evolution groups to disentangle the effects of normal evo- 
lution abundance changes (Section I4.2j) . 

4.1. Mean Abundances of the Stellar Samples 

The first question we wish to address is whether the 
abundance distribution of light elements differs between 
the slow and rapid rotator samples. We plot the Li abun- 
dance as a function of projected rotational velocity in 
Figure [TO] This plot gives us the first evidence from 
our own data that there is a relationship b etwee n rota- 
tional velocity and A(Li), as IDrake et al.l (|2002f ) previ- 
ously suggested. The average ^4(Li) and standard errors 
in the mean for the entire slow rotator and rapid rota- 
tors samples are —0.18 ± 0.08 and +0.81 ± 0.27, respec- 
tively, and these mean values (without the errors) are 
indicated on the plot. The average abundances include 
the upper limit measurements. If the SWPs are excluded 
from the slow rotator sample, the mean abundance in- 
creases to — 0.14± 0.09. Five of the fifteen rapid rotators 
have j4(Li) exceeding the maximum ^4(Li) of all of the 
slow rotators. A two-sided Kolmogorov Smirnov (K-S) 
test reveals that the probability that the lithium abun- 
dances of the slow and rapid rotators are from the same 
parent population is only 0.6%. We can subdivide the 
rapid rotators into the two RV groups defined in Section 
13.3.21 the CB group, which contains the stars that may 
be synchronously rotating with binary companions, and 
the SLP group, which require some other explanation 
(such as planet accretion) for the rapid rotation. The 
CB group has a mean A(Li) of +0.41 ± 0.37 dex, which 
is much closer to the mean A(Li) of the slow rotators 
than the A(Li) of the entire ensemble of rapid rotators 
is to the slow rotators. The CB group is small enough 
that it has a comparable probability of being drawn from 
either the slow rotator sample (48%) or the RV stable, 
rapid rotator sample (53%). The SLP group, in contrast, 
has a mean ^4(Li) of +1.06 + 0.36 dex, which is even more 
Li enhanced than the entire rapid rotator ensemble. The 
K-S probability is 0.5% when comparing the distribution 
of ^4(Li) of the SLP group to the slow rotators. 

We make a similar comparison of 12 C/ 13 C in our stel- 
lar samples by plotting 12 C/ 13 C against vsmi in Figure 
[TT1 The average 12 C/ 13 C of the slow and rapid rotators 
subsamples are overplotted. The rapid rotators again 
seem to show the enhancement expected from planet ac- 
cretion; they have 12 C/ 13 C = 19.0 + 1.9 compared to the 

slow rotators that have 12 C/ 13 C = 17.0 ± 0.9. By again 
subdividing the rapid rotators, we find that the nine stars 
in the SLP group have 12 C/ 13 C = 17.3 ± 2.2. Only four 
of the CB stars have measurable 12 C/ 13 C, and the av- 
erage of those stars are 23.0 ± 3.4. However, because of 
the larger relatively uncertainties in the 12 C/ 13 C mea- 
surements, K-S tests of the 12 C/ 13 C distributions show 
no statistically significant differences between the slow 
rotators, SWPs, RRs, and subsamples of RRs. 

4.2. Common Evolution Groups 

Both the strength and weakness in the above compar- 
ison is that it averages over the normal variations in the 
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Figure 10. Plot of A(Li) as a function of vs'mi for the slower 
rotators with and without known planets (triangles and squares, 
respectively), rapid rotators that are single or in long-period bi- 
naries (SLP group, circles), and rapid rotators suspected to be in 
close binaries (CB group, x). The horizontal dashed lines show 
the average abundances of the slow (i>sini< 8kms _1 ) and rapid 
(vsini> 8kms -1 ) rotators. Upper limits are indicated with down- 
ward facing arrows. The sizes of the typical errors are shown in 
the upper left corner. 
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Figure 11. Plot of 12 C/ 13 C as a function of jjsini for the slow 
and rapid rotators, using the same symbols as in Figure [TOl The 
dashed line indicates the mean 12 C/ 13 C of the slow and rapid 
rotators. Lower limits are indicated with upward facing arrows. A 
typical error bar is shown on one of the slow rotators. 

abundances to probe how rotational velocity is poten- 
tially correlated to the abundances. However, there are 
well-known changes that are expected in these elements 
as the stars evolve through the red giant phase. Our slow 
and rapid rotators are not distributed equally across the 
various stages of RGB evolution; therefore, it is worth re- 
peating the comparison of -A(Li) and 12 C/ 13 C with v smi 
in subgroups of similarly evolved stars. 

We define evolution groups by comparing our data 
to model ev o lution tracks. In Figure Q21 we plot 
iGirardi et al.l (|2000f ) solar metallicity stellar evolution 
tracks on a modified HR diagram (replacing luminos- 
ity with logg) for the mass range of 0.8-3.5 M . The 
color scale of the evolution tracks indicate the relative 
depth of the convection zone ranging from the shal- 



lower depth the convection zones have on the subgiant 
branch to the full FDU depth. Major stages of stellar 
evolution — the first ascent RGB phase, the luminosity 
bump, and the post RGB-tip evolution — shift to cooler 
temperatures and larger surface gravities at super-solar 
metallicities (and vice versa for sub-solar metallicities) . 
The rapid rotator sample is small enough that grouping 
stars first by similar metallicity and then by evolution 
stages would generate subsamples too small to make use- 
ful comparisons. Instead, to compare our stars to these 
solar-metallicity tracks, we offset the observed T e g and 
log<? of our program stars by metallicity-dependent fac- 
tors to create T e ' ff = T cff + [Fe/H] x SOOKdex" 1 and 
log g' = log g — 0.25 x [Fe/H]. These "primed" stel- 
lar parameters represent the temperatures and gravities 
the target stars would have if they were so lar-metallicity 
stars. The offsets were found by comparing lGirardi et al.l 
(2000) evolution tracks of Z = 0.004, 0.008, 0.019, and 
0.030. Stellar evolution groups are defined by making 
cuts perpendicular to the general direction that stars 
evolve through the T e ' ff -logf/ plane, preferably in loca- 
tions where there were gaps in the distribution of our 
program stars. These cuts were done by eye, and are in- 
dicated by the dashed lines in Figure [T^l One additional 
cut, indicated by the dotted line, was made to separate 
stars that are at the luminosity bump from those stars 
have not yet reached this important stage. None of the 
rapid rotators fall into the luminosity bump class. 

The groups are labeled according to the normal stellar 
evolutionary processes affecting the j4(Li) and 12 C/ 13 C 
of solar-mass stars at that stage. These groups (from 
bottom left to upper right in Figure fT2]) are as follows. 

Dilution in Progress.: Stars are still undergoing Li di- 
lution. These stars should have the largest A(Li) 
and 12 C/ 13 C in normal stellar evolution. 



Significant Dilution.: 

maximum depth; 
tially complete. 
12 C/ 13 C is falling. 



The convection zone is near its 
therefore, Li dilution is essen- 
13 C is being dredged-up, and 



Max CZ Depth.: The convection zone (CZ) has 
reached its maximum depth and is receding. A (Li) 
and 12 C/ 13 C are at standard dilution values. The 
H- burning shell is still interior to the maximum CZ 
depth. 

Bump.: The H-burning shell reaches the maximum CZ 
depth, and non-convective mixing can occur. 7 Be 
created by the Cameron Fowler chain may be car- 
ried into the convection zone where it decays to 
7 Li, which is not destroyed. There is a potential 
for large A(Li) at this stage, perhaps even exceed- 
ing that seen in the "Dilution in Progress" group. 
12 C/ 13 C remains constant at the standard value. 

Early Post-Bump.: Stars have evolved just beyond 
bump. Non-convective mixing may be reducing 
both A(Li) and 12 C/ 13 C. Li can be either rather 
high (if Li was briefly regenerated) or very low (if 
no Li was regenerated). 12 C/ 13 C may be dropping. 

Late Post-Bump.: Low to very low values of A (Li) and 
12 C/ 13 C depending on the amount of extra mixing. 
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Figure 12. Illustration of our selection of stellar evolution groups. The "primed" T e ff and logg are the observed values plus a metallicity- 
depen dent shift applied to align them to solar metallicity evolution tracks (see the text for details). The thin solid lines are lGirardi et aLl 
(2000) solar metallicity evolution tracks for masses between 0.8 Mq (bottom track) and 3.0 Mq (top track). Evolution progress along 
each track from lower left to upper right. The color scale of the models illustrates the relative depth of the convection zone (CZ) from 
shallow (light gray or red) to deep (dark gray or blue). The heavy solid line marks the stage at which the maximum CZ depth is achieved. 
Our stellar sample is plotted with the same symbols identified in Figure ITol The dashed lines were chosen by eye to group stars of similar 
evolutionary stages. The dotted line is used to isolate the group of stars that fall on the luminosity bump. The bump can be identified as 
where lower mass tracks briefly double back in the direction of evolution. 



In Figure fT3l we plot .A (Li) versus 12 C/ 13 C for these six 
stellar groups along with two models of how these prop- 
erties should change in stars of different masses. The 
first model star is 3M Q , for which only sta ndard dilu- 
tion is expected to affect the abundances ( Ibenj Il967t 
lLambert et al.l fl980( ). In this model, the A(Li) at the 
end of the MS evolution is ~ 3 dex. The second 
model is for a 0.85 M Q star for which standard dilu- 
tion and canonical extra-mixing a re expected to occur 
(|Denissenkov fc VandenBerd 120031) . The A(Li) at the 
end of the MS for this model is ~ 2.4 dex. We expect 
the stars in the first three evolution groups to fall be- 
tween the two model lines in the upper right portion of 
each panel. For the remaining three evolution groups, we 
expect any higher mass stars in our sample to remain in 
the upper right region of the plot (near the dotted line), 
while lower mass stars should move down and left along 
the low mass model. Our data, however, do not follow 
these expectations. In all of the subgroups, there are no 
clear correlations between A(Li) and 12 C/ 13 C (though, 
one should remember that the errors in 12 C/ 13 C are of- 
ten large: typically ±5 for 12 C/ 13 C> 20). Almost all of 
the stars appear to show substantial Li depletion (their 
Li abundances fall under the 3 M Q model), and at the 
coolest temperatures ("Late post-bump" stage), some of 
the stars have Li below the extra-mixing model. More 
than half of the stars at each stellar evolution stage have 
12 C/ 13 C that is lower than the standard model (i.e, lie 
to the left of the 3 M© model) , but none of the stars 
have a 12 C/ 13 C (within the measurement errors) that is 
lower than that predicted by the extra-mixing model. If 



real, the unexpectedly low 12 C/ 13 C at early RGB stages 
could indicate a variation in the initial 12 C/ 13 C of the 
stars. Similarly, because the variation in A (Li) at each 
evolutionary stage exceeds the measurement errors, we 
are likely seeing intrinsic scatter in A(Li) that makes a 
simple comparison to only two models difficult. A(Li) 
that is lower than model predictions can be explained by 
the variations of slow Li-depletion processes on the MS. 
Both models use MS A (Li) abundances that are at least 
an order of magnitude larger than the present-day solar 
values. On the other hand, the two Li-rich rapid rotators 
in the "Significant Dilution" group exceed the predicted 
values even if no MS Li-depletion has occurred. 

To compare the effect of rotation on A (Li) more ex- 
plicitly, we plot A(Li) against vsmi for each of the stel- 
lar evolution bins in Figure [TJJ The average A (Li) of 
the slow and rapid rotators in each bin are indicated 
with lines. This figure is meant to address the question 
of whether the larger mean A(Li) of the rapid rotators 
(discovered in Figure [T0| was merely result of the rapid 
rotators tending to be less evolved on average compared 
to the slow rotators. It is noteworthy, therefore, that 
in all of the evolution bins that contain more than one 
rapid rotator, the mean A(Li) of the rapid rotators ex- 
ceeds that of the slow rotators. Furthermore, the most 
Li-rich star in the bin is a rapid rotator. Curiously, in 
both evolution bins that contain fewer than two rapid 
rotators, there are a distinct group of slow rotators with 
A(Li) larger than the rest. 

5. DISCUSSION 
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Figure 13. A(Li) plotted against 12 C/ 13 C for the slower rotators (triangles), giants with known planets (squares), the SLP group of rapid 
rotators (circles), and the CB group (x). Each panel shows a stellar evolution group, progressing from less evolved to more evolved (left to 
right, then top to bottom). Upper limits in A(Li) and lower limits in 12 C/ 13 C are indicated with downward and rightward facing arrows, 
respectively. The lines show the expected relation ship of these two properties for a 3 Mq star (dotted), for which only standard dilution 
is expected (e.g., Ibcn 1967; Lambert ct al. 1980), a nd for a low- mass star (0.85 Mq, dashed line), for which both standard dilution and 
canonical extra- mixing are expected (Dcnissenkov & VandcnBcrg 2003). 



In the previous section, we found two important re- 
sults. First, the rapid rotators tend to be more Li-rich 
than the slow rotators regardless of evolutionary stage. 
Second, there are at least two stars with A(Li) exceeding 
even the most conservative estimates of Li dilution (the 
standard dilution models) as well as a handful of stars 
that appear to be Li-rich when comparing their abun- 
dances to the other stars of similar evolutionary stage, 
but whose abundances compared to the models in Figure 
1131 are not atypical. Here, we consider various scenarios 
for understanding these two results. 



while stars with more modest Li enrichments will show 
abnormally low 12 C/ 13 C. On average, therefore, 12 C/ 13 C 
should be lower in stars that have undergone Li regener- 
ation. Considering these two observational signatures of 
Li regeneration, we can conclude that our Li-rich rapid 
rotators do not fit this picture for the following reasons. 
First, on average the rapid rotators do not exhibit lower 
12 C/ 13 C than the slow rotators. Second, if our stellar 
evolution groups are defined correctly, then the two most 
Li-rich stars are at pre-bump stages and should not have 
yet experienced Li regeneration. 



5.1. Post-bump Lithium Regeneration 

The prevailing explanation for Li-rich stars is inter- 
nal Li regeneration via the Cameron- Fowler chain, which 
is thought to begin at the luminosity bump when the 
convective envelope can be connected to the H-burning 
shell via non-convective (generally diffusive) mixing pro- 
cesses. At the luminosity bump, the hydrogen burning 
shell erases the chemical discontinuity, or "//-barrier," 
that inhibits mixing. Considerable work has been put 
into trying to understand the physical causes of this 
"cool-bottom processing" (CBP) mechanism — whether 
the mixing results from, e.g., a thermohaline instability 
or buoyant magnetic flux tubes. Regardless of the under- 
lying physics, one observational consequence of this ex- 
planation is that the star must have evolved to the lumi- 
nosity bump, and indeed, iCharbonnel fc Bala chandran 
(2000) found most Li-rich giants to be associated with 
this evolutionary stage. Another observational signature 
of this processing is the further reduction of 12 C/ 13 C 
that occurs after the Li regeneration ceases. Stars with 
the largest A (Li) may appear to have normal 12 C/ 13 C, 



5.2. Rotation-induced Mixing 

An alternative Li-regeneration solution is one in which 
the enhanced rotation of the rapid rotators drives the 
non-convective mixing that connects the envelope to the 
Li-burning interior. The assumption is that the rapid 
rotation creates favorable mixing currents that can de- 
posit the newly regenerated Li (or the parent isotope 
7 Be) into the stellar envelope. In this scenario, the star's 
evolutionary proximity to the luminosity bump is irrel- 
evant because the enhanced rotation is itself responsi- 
ble for the non-canonical mixin g. Rotation-induced mix- 
ing ha s be en explored by both iSackmann fc Boothrovdl 
(1999) and lDenissenkov'&: Herwid (|2004f ) to explain Li- 
rich giants, and both models are also capable of lowering 
12 C/ 13 C from the standard model predictions. However, 
both groups were also exploring rotation-induced mix- 
ing as a possible CBP mechanism. In other words, their 
models required the erasure of the chemical discontinu- 
ity at the luminosity bump, whereas a successful model 
for explaining our pre-bump Li-rich rapid rotators is one 
which does not require the erasure of the mean molecular 
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Figure 14. Same data as Figure 1101 but plotted for the stellar eve 
right, then top to bottom). The symbols have the same meaning as i 
the slow rotators (left line) and rapid rotators (right line). 

weight discontinuity to function. IChaname et al.l (|2005l ) 
were the first to relax the assumption that the /x-barrier 
inhibits all extra mixing to see whether rotational mixing 
would occur at other evolutionary phases. They explored 
the evolution of different rotational profiles and found 
that rotational mixing does indeed become active beyond 
the luminosity bump but only if the CZ was differentially 
rotating. In other words, in the case of solid body rota- 
tion in the stellar envelope, a star would not experience 
extra mixing even after the /u-barrier was erased. Addi- 
tionally, they found that the FDU mixing of 12 C/ 13 C oc- 
curred earlier and more gradually in their rotating mod- 
els compared to standard models (though both models 
resulted in the same post-FDU value). On the other 
hand, to reproduce the abundance observations, the ini- 
tial rotation rates of th eir models must e xceed the values 
actually observed. The IChaname et al.l (|2005[ ) paper did 
not explore the explicit effect of rotation on the Li abun- 
dance, in particular whether their models allow the spe- 
cial mixing cases capable of creating at least short-lived 
periods of enhance Li (as opposed to further destruction 
of Li that generally occurs with extra mixing). Further 
work is needed to see whether such conditions are possi- 
ble and whether those conditions require the removal of 
the "/x-barrier" or not. 

5.3. Helium Flash 

Recently, iKumar et al.l ()2011l . hereafter, Kll) con- 
ducted a large survey for Li-rich stars in a sample of 
2,000 red giants and found 15 previously unidentified Li- 
rich stars. Like our study, they find some Li-rich giants 
that are too warm to be associated with the luminos- 
ity bump. Kll suggest that these stars are likely as- 
sociated with the red clump — the core helium-burning 
stage of relatively metal-rich stars. They suggest two 
scenarios for explaining the presence of large ^4(Li) in 



ition g roup s progressing from less evolved to more evolved (left to 
Figure [131 The dashed lines in each panel show the mean A(Li) for 

post-RGB stars. The first explanation is simply the 
survival of Li that was regenerated at the luminosity 
bump; however, this scenario is contradicted by the rel- 
atively few Li-rich stars between the bump and RGB 
tip. Kll also hypothesize the possibility that Li could 
be regenerated during the He flash, although this sug- 
gestion is based simply on the observed coincidence of 
their Li-rich giants near the red clump more so than on 
any physical mechanism. The idea is that some frac- 
tion of the 3 He remaining in the RGB tip star could 
be converted to Li through the Cameron-Fowler process, 
whereby 3 He( 4 H,7) 7 Be(e~, v) 7 Li. They suggest that this 
model would only operate in stars within the narrow 
mass range of 1.5Mq < < 2.25M Q , where the star 
has a low enough mass to experience a He-flash but has 
a high enough mass for some 3 He to survive. As a point 
of interest, we believe that Kll put a conservativ e lower 
limit on the mass range. Eggl eton et al.l ([2008) mod- 
eled the destruction of 3 He and found that as much as 
~ 95% of the 3 He created during the MS evolution will 
be destroyed by the time the star reaches the tip of the 
RGB. For the 1.5 Mq model, the destruction fraction 
is only 75-83%. Nevertheless, the a mount of 3 H e cre- 
ated in low mass stars is substantial. Ilbenl (|1967f ) esti- 
mates ~ 10 -3 Mq, which corresponds to ^4( 3 He) ~ 8.6. 
Assuming 95% of the original 3 He was destroyed, only 
10~ 4 of the remaining 3 He must end up as 7 Li to yield 
A(Li)— 3.3 dex (the abundance of our most Li-rich star). 
In the 1.5 Mq case, the fraction of 3 He remaining is 
~ 20 %, but the original 3 He budget is smaller, so that 
the reservoir of 3 He at the RGB tip only about 2.5 times 
larger in the 1.5 Mq case than in the 1 Mq case. Kll do 
not report v sini, so we cannot comment on any relation- 
ship of rapid rotation to -A(Li) in their results. However, 
Kll do provide 12 C/ 13 C, and in contrast to our stars, 
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their warm Li-rich giants tend to have 12 C/ 13 C< 16, in- 
dicating that some sort of extra mixing process has oc- 
curred in their Li-rich stars that does not appear to have 
occurred in the Li-rich stars studied here. 

5.4. Planet Accretion 

One of the main drawbacks to the Li-regeneration mod- 
els is that they fail to explain the excess angular momen- 
tum of the rapid rotators. The accretion of a planet is 
one means by which a red giant star can acquire sufficient 
angular momentum to become a rapid rotator. We can 
test whether our abundance results are consistent with 
the planet accretion paradigm by estimating the masses 
and chemical compositions that accreted planets would 
have had to account for the mean abundance differences 
between the slow and rapid rotators. The expected stel- 
lar abundances of Li after planet accretion — A(Li) new — 
are given by 



A(Li), 



log(g o 10^ Li ^ + 10^*) - log(l + ?.), (2) 



where q e is the ratio of the planet mass (M p ) to the mass 
in the stellar convective envelope (M env ) and A(Li) p and 
.A(Li)* are the initial abundances of Li in the planet and 
star, respectively. (See the Appendix for the deriva- 
tion of this equation.) For the present argument, we 
assume that the slow rotators are representative of the 
initial stellar abundances (i.e., A(Li)* = ^4(Li) slow = 
— 0.18dex) while the RV stable rapid rotators represent 
the post-planet-accretion abundances (i.e., v4(Li) ncw = 
1.06 dex). If we assume that the meteo ritic Li abundance 
of ou r solar system (A(Li) ~ 3.3 dex, lLodders fc Fegleyl 
fl998l) is a good representation of the Li abundances of 
hypothetically accreted planets (i.e, A(Li) p = 3.3), then 
we find from Equation @ that q c — 5.4 x 10~ 3 . This q e 
corresponds to ~ 6Mj up assuming M cnv = 1 M Q , which 
is a reasonably-sized planet accreted into a reasonably- 
sized red giant convection zone. 

Similarly, we can test our non-detection of a 12 C/ 13 C 
increase using our estimated q c . The carbon ratio ex- 
pected after the accretion of a planet, ( 12 C/ 13 C) nC w, is 
given by (see the Appendix) 



( i2 C/ i3 C) 



10^(c)p^ + io^(c), 
i+i-p 

1QA(0) P ^_ + 10 A(C) 4 



1+r- 



1+T-* 



(3) 



where r p = ( 12 C/ 13 C) P and r* = ( 12 C/ 13 C%. For a solar 
metallicity red giant star, we take A(C)* = 8.26, which is 
the s olar-metallicity (A(C)q = 8.39 dex; iGrevesse et al.l 
120071) adjusted for some po st-MS processing of C a nd N 
such that [C/Fe] = -0.13 (|Gustafsson et al.|[200l . As 
with the lithium example, we use the average slow ro- 
tator 12 C/ 13 C to represent pre-planet accretion so that 
= 17.0. Using Jupiter as an analog, we expect 
the assimilated planet to have ha d 3x the solar car- 
bon abundance (|Wong et al.l 120041) . Thus, A{C) P = 
A(C) Q + log(3) = 8.87 dex. We adopt r D = 89, which is 
a sta ndard value in the solar system (jLodders fc Feglevl 
1998). Under these assumptions and adopting the q e de- 
rived from the Li enhancement, we calculate that a sam- 
ple of planet accreting stars (such as the rapid rotators) 
is expected to ha ve 12 C/ 13 C =17.3. Saturn is even more 
carbon-enriched (jMousis et a l. 2009) with up to 10 x the 



solar value. Using Saturn's composition as the analog 
for accreted planets, the expected carbon abundance in- 
creases to 12 C/ 13 C = 18.0. The average 12 C/ 13 C of the 
RV stable rapid rotators (17.3) is similar to the value ex- 
pected from accreting a Jupiter analog, supporting the 
plausibility of planet accretion by these stars. However, 
given the relative sizes of the signal (the 12 C/ 13 C in- 
crease) compared to our uncertainties, we caution that 
we cannot conclusively assert that we are measuring a 
true 12 C/ 13 C difference. 

We can further test the plausibility of the planet ac- 
cretion hypothesis by assessing whether both the en- 
hanced angular momentum an d Li enrichment are con- 
sistent. Carl berg et al.l (|2009l ) introduced an equation 
(their Equation (1)) to relate an observed stellar usini to 
the properties of the star and the planet it accreted, i.e., 
v sin i = 8 [M p sin i y/GM+a^l - e 2 )] /M env R* , where M p 
is the planet mass, a p is the planet's initial orbital sep- 
aration, and -Me nv is the mass in the stellar convection 
envelope. The terms in brackets describe the initial angu- 
lar momentum of a planet orbiting the star. To estimate 
the angular momentum gained by the star, both the ap- 
proximate stellar envelope mass (M env ) and radius must 
be known (i?*). Because the latter can change by two 
orders of magnitude during RGB evolution, it is helpful 
to perform the analysis on individual stars as opposed 
to using the average stellar properties we have used thus 
far. To this end, we select a small group of Li-rich rapid 
rotators that represent the best candidates in our sam- 
ple of stars that have accreted planets. In addition to the 
enhanced rotation and enriched Li, these are stars that 
have not yet evolved to the luminosity bump and are 
not suspected to be in close binary systems. Three stars 
meet all of these requirements and are listed in Table [12j 
The first two stars in that table are the two most Li- 
rich stars in the "Significant Dilution" group, as seen in 
FigureQH The third star, Tyc3340-01195-1, is the most 
Li-rich star in the "Dilution in Progress" group and has 
a long-period binary companion. 

Matching the stars' effective temperatures, surface 
gravities, and metallicities to grids of either stellar evo- 
lution tracks or isochrones can be used to estimate the 
stellar masses and radii. We opt to use isochrones in this 
analysis because the isochrones of iMarigo et al.l (2008) 
can be interpolated onto a finer grid in Z than what is 
available for the evolutionary tracks. We downloaded 13 ! 
isochrones with Z = 0.0002^ 0.001, 0.004, 0.008, 0.012, 
0.016, 0.020, 0.024, and 0.03 (Fe/H = -2.05, -1.33, 
-0.72, -0.41, -0.23, -0.10, 0.00, +0.08, and +0.18), 
and logt ranging from 8.5 to 10.1 in increments of 0.1 dex 
(where t is the age of the stellar population in years). 
The stellar masses and radii of the stars are estimated 
in the following manner. First, we identify the isochrone 
metallicity that most closely matches the observed stellar 
metallicity. For each age at constant Z, we find where the 
isochrone intersects the stellar T e g. If there are multiple 
intersections, which is common for red giant evolutionary 
stages, we use the spectroscopically derived log g to se- 
lect the best intersection. The isochrones are computed 
on finite grids of T c ff points, so we identify the isochrone 
points that straddle the stellar T e s and interpolate the 

13 http://stev.oapd.inaf.it/cmd accessed on 2011 May 17. 
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Table 12 

Best Case Candidates for Planet Accretion 



Star 




v sin i 


A(Li) LT E 


^(Li)nlte 


12 C/ 13 C 


[Fc/H] 




(K) 


(kins -1 ) 


(dcx) 


(dex) 




(dcx) 


G0928+73.2600 


4770 


8.4 


+3.62 


+3.30 


28 ± 8 


-0.26 


Tyc0647-00254-l 


4825 


10.4 


+ 1.92 


+2.06 


20 ± 3 


-0.01 


Tyc3340-01195-1 


5040 


8.4 


+ 1.21 


+ 1.32 


25 ± 5 


-0.18 



3.0 



1.0 




1 2 3 4 5 6 7 

M (M s ) [Lit.] 

Figure 15. Comparison of the isochrone-derived masses of the 
giant SWPs (j/-axis) with the mass measurements reported in the 
Exoplanet Encyclopedia (x-axis). A unity-slope line is plotted for 
reference. 

masses and log g associated with those two points to the 
T G ff we are interested in. In other words, for each star 
we reduce the isochrones to a grid of age, mass, and 
log<? for the stellar temperature. We then find adjacent 
grid points of isochrone-derived \ogg that straddle the 
stellar log g, and we interpolate the ages and masses as- 
sociated with those points to estimate the age and mass 
of the stellar log g. This last step generally finds be- 
tween one and three unique age/mass solutions, and the 
final mass and age estimates of our program stars average 
over these unique solutions. Once we estimate the stellar 
mass, it is trivial to calculate an estimate in radius. Us- 
ing the measured surface gravity, the radius is given by 
R* = \J GAI*/ g, where G is the gravitational constant. 
Although we acknowledge the large uncertainty inherent 
in this analysis, especially with the overlap of the RGB 
and horizontal branch at the stars' T e ff and log g, we find 
that our isochrone fitting works rather well. The sample 
of giant SWPs has independent measurements of mass 
with which we can compare our results. In Figure H~5l we 
plot the comparison of our mass measurements to litera- 
ture mass measurements and find that our masses agree 
with literature values within the quoted uncertainties, 
with the exception of Pollux, which has Mut. — 1.9M 
and M isoch = 2.2 ± 0.2 M Q . 

The isochrones do not provide information on the in- 
ternal structure of the stars. Therefore, we use the 
IGirardi et al.l ([20000 stellar evolution models that most 
closely match the mass and metallicity of our stars to 
estimate the fraction of the mass in the convective en- 
velope. The least evolved star in Table [T^] is Tyc3340- 
01195-1, which has ~ 32% of the total mass in the stellar 
envelope. The other two stars have envelope mass frac- 
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Figure 16. Minimum "projected mass" (M p sini) needed to ac- 
count for the observed v sin i of the best case planet accretion can- 
didates as a function of the expected maximum initial semimajor 
axis (ao) of the accreted planet. The true mass of any accreted 
planets will be larger both for smaller inclination angles (more 
pole-on angles) and for ao < max.(ao). For reference, the plan- 
etary masses and orbital separations of known extrasolar planets 
orbiting MS stars are shown with crosses. 



tions of - 80%. 

Obviously, the unknown inclination angle (i) limits 
us to estimating the minimum planet mass (M p sini). 
There is also a degeneracy between the planetary mass 
and orbital separation. A massive planet at a small sepa- 
ration can create the same degree of stellar rotation as a 
less massive planet at an initi ally larger separation. For- 
tunately, the tidal modeling of lCarlberg et al.l (2009) pro- 
vides a means of breaking this degeneracy in the limiting 
case of the maximum initial orbital separations of planets 
that could have been accreted (a max ) — planets that are 
too distant from their stars will n ot be accreted. In the 
exoplanetary systems modeled in iCarlberg et al.l (|2009t 
99 systems with a total of 115 planets), we find that the 
ratio of the planets' initial semimajor axis (ao) to the 
stellar radius at the time of planet accretion, i?*(t a cc), is 
typically 4.2. Therefore, we can estimate a max ~ 4.2 i?*. 
Substituting this maximum orbital separation into the 
angular momentum equation yields the minimum planet 
mass capable of producing the observed rotational ve- 
locities in the rapid rotators. This derived mass is both 
a minimum M p because of the unknown inclination and 
a minimum M p sin i because more massive planets with 
smaller initial a p are also capable of producing the ob- 
served rotation. The results of this calculation for our 
selected Li-rich rapid rotators are presented in Figure 
1161 which shows the minimum M p sin i required to spin 
up the star as a function of a max ~ 4.2i?*. The least 
evolved star could have accreted planets that originally 
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orbited within 0.1 AU and requires a Jupiter- mass planet 
to account for the angular momentum. The other two red 
giants could have accreted planets within 0.2 AU, and the 
planets must have been at least 4 Mj up to explain their 
enhanced rotation. This minimum planet mass is slightly 
less than the average of 6 Mj up planet needed to account 
for the Li-enrichment in Section 14.11 For comparison, 
we also plot in Figure [16] the planetar y m asses and or- 
bital separations of known exoplanetsQ Our angular 
momentum analysis implies that our rapid rotators can 
be explained by accreting planets to the left and above 
the rapid rotators' positions on the plot, and we find such 
planet do exist among the known exoplanetary systems. 
However, it is also worth noting that there is a dearth 
of Jupiter mass planets (and larger) between ~ 0.08 and 
~ 0.7 AU, the region where stars like those listed in Table 
[T2l would be actively clearing out their planets by tidal 
decay. 

5.5. A Few Noteworthy Stars 

G0928+73.2600 This star was analyzed separately 
bv lCarlberg et al.l (|2010bl ). with attention brought to its 
apparent pre-bump evolutionary stage but high enough 
A(Li) to require accreting a Li-enriched object. We high- 
light this star again to contrast it to another star in the 
study that has almost identical stellar parameters, but 
very different A(Li) and 12 C/ 13 C. Using the stellar pa- 
rameters from the 2007 observations, G0928-I-43.2600 has 
T cff = 4900 K, log5= 2.70, and [Fe/H]= -0.26, and the 
slow rotator G1200+67.3882 has T cff = 4900 K, log g= 
2.70, and [Fe/H]= —0.24. Aside from rotation speeds 
(8.4 kms -1 and 1.5 kms -1 ), these stars differ in their 
abundances. The rapid rotator has A(Li)= 3.30 and 
12 C/ 13 C= 28, while the slow rotator has A(Li)= +0.39 
and 12 C/ 13 C= 8.6. Such a low 12 C/ 13 C is suggestive 
of more than the usual mixing, which implies that this 
star could have regenerated Li internally and is now in 
the process of destroying that regenerated Li and reduc- 
ing 12 C/ 13 C below the normal values. If this scenario 
is indeed true, it implies that G0928+73.2600 is at the 
slightly earlier stage where newly regenerated Li has not 
yet been destroyed and so the 12 C/ 13 C is still at the 
standard value. 

The problem with using Li regeneration to explain the 
abundance differences of these otherwise similar stars 
is that the stars are both too hot to have evolved 
through the bump phase (T bump ~ 4600-4700 K for 
[Fe/H]=-0.25). This difficulty could be resolved if 
our choice of stellar evolution tracks predicts too cool 
temperatures for the luminosity bump or if our de- 
rived spectroscopic temperatures are systematically too 
hot. To test the first scenario, we compare stel- 
lar evolution tracks from t wo in dependent sources. 
ICharbonnel fc Balachandran (2000) plot solar metallic- 
ity stellar evolution tracks in their Figure 1. A visual in- 
spection of that plot shows the luminosity bump at tem- 
peratures between 4450 and 4500 K, which is ~ 100 K 
cooler than the luminosity bump of the solar-metallicity 
models we plot in Figure [121 As a second test, we com- 
puted the evolution of a 2 M Q , Z = 0.011 star us- 

14 D ata come from the Exoplanet Orbit Database IWright et al.l 
1)20111 ). accessed on 2011 December 18. 



ing the MESA code (jPaxton et al.l 1201 llPI . The lumi- 
nosity bump of that model extends to a temperature 
of 4700 K, comparable to the models shown in Figure 
1121 Next, we compile two different photometric temper- 
atures to test for a systematic offset in our spectroscopic 
derivation. The firs t set of temperatures comes from 
iCarlberg et ail pOll and were derived using the stars' 
Washingt on M—Ti colors ( which are co nverted to Cousin 
V - I , iMaiewski et~aTl 120001 ) and iHoudashelt et al.l 
(2000) color-temperature relations. This calculation 
yields T p hoti = 4773 K for the rapid rotator and Tphoti — 
4875 K for the slow rotator. Both stars also have Ty- 
cho designations. G0928+73.2600 is Tyc4382-00780-l, 
and G1200+67.3882 is Tyc4160-00999- l. We converted 
their observed Tycho B — V magnitudes (|Hpg et al.ll2000D 
to Johnson B — V = 0.85( . B — U )Tvnho, dereddened the 
colors using iSchlegel et all (fl99l maps, and again used 
IHoudashelt et all ([20001 ) color-temperature relations to 
find T p hot2 = 4755 K for the rapid rotator and T p hot2 = 
5002 K for the slow rotator. Thus, the rapid rotator 
may be slightly cooler and more evolved than our spec- 
troscopic analysis suggested, while the slow rotator may 
be slightly warmer (and less evolved). However, the pho- 
tometric temperatures are still warmer than the luminos- 
ity bump temperatures. Therefore, these two stars may 
represent examples of the He-flash lithium regeneration 
hypothesized by Kll. 

Tyc3340-01195-1. This star is a long period binary 
star, and Cll argued that the relatively large sepa- 
ration of the stellar components in the former system 
(a* sin i ~ 425i?o ~ 2 AU) made tidal synchronization 
an unlikely explanation for the primary star's enhanced 
rotational velocity (v sini=8.4kms -1 ). The presence of 
a stellar companion raises questions about the stability 
of planets in the system. However, planets can have sta- 
ble orbits around the primary star interior to the stellar 
companion if the orbits are smal l enough. In a s t udy o f 
such "S-type" planetary orbits, iRabl fc Dvorak! (|1988f ) 
found a quadratic relationship relating the largest Low- 
est Critical Orbit (LCO — orbits larger than this may 
be unstable, while smaller orbits are stable) to the bi- 
nary system's separation (a*) and eccentricity (e) to be 
LCO = 0.262a* -0.254a*e-0.060a*e 2 . For the Tyc3340- 
01195-1 system, a planet orbiting with a p < 0.4 AU 
would be stable around the primary. In Figure[16]we find 
that the maximum initial separation of planets that could 
have been accreted by Tyc3340-01195-1 is only - 0.1 AU; 
therefore, the orbits of any accreted planets would have 
been dynamically stable in the binary star system during 
the primary star's MS lifetime. 

6. SUMMARY 

We have studied the global abundance patterns of 
A(Li) and 12 C/ 13 C in a homogeneously selected sample 
of slow and rapid rotators to see whether the rotation is 
correlated with the replenishment of elements destroyed 
during the stellar evolution, as expected from planet ac- 
cretion. Our final sample contains 71 slow rotator stars 
(58 chosen purely for their slow rotation, 10 selected be- 
cause they are known to host planets, and the three M08 

5 MESA version 3661, using the "inlist" file provided in the 
"lM_pre_ms_to_wd" test suite as a template. We changed only the 
mass and metallicity of that file. 



2G 



Carlberg et al. 



stars) and 15 rapid rotators. From an analysis of line- 
of-sight considerations, we expect that 4.5% of our slow 
rotators have true rotational velocities qualifying as rapid 
rotation. 

Overall, the rapid rotators show an enhancement of 
A(Li) over the slow rotators by 0.99 dex. When select- 
ing the subset of rapid rotators that are RV stable, the 
average enhancement of ^4(Li) increases to 1.24 dex over 
the slow rotators. These Li enhancements are consistent 
with the accretion of a ~ 6 Jupiter masses of planetary 
material with a Li abundance similar to meteoritic abun- 
dances (thought to be relic of the solar nebular abun- 
dance). Consistent with this explanation and our rela- 
tively large 12 C/ 13 C uncertainties, we find no statistically 
significant difference between the 12 C/ 13 C measured in 
the rapid and slow rotator samples. A more massive or 
more carbon-enriched object would have to be accreted 
to be measurable at the level of our 12 C/ 13 C uncertain- 
ties. However, we do not measure smaller 12 C/ 13 C in 
our Li-rich stars or rapid rotators, a signature that would 
need to be seen to prove the replenishment of Li through 
nuclear processing. We also compared our stellar sam- 
ple to evolutionary tracks to ascertain how the relative 
Li and carbon abundances varied within subsets of sim- 
ilar evolutionary stage. This comparison was necessary 
because the rapid rotators and slow rotators are not dis- 
tributed evenly across the RGB. We found that in all 
groups that contained more than one rapid rotator, the 
most Li-rich star in the group was a rapid rotator. In 
other words, the result that rapid rotators are more Li- 
rich than the slow rotators persists even within groups 
of similarly-evolved stars. 

These main conclusions were drawn by comparing the 
global properties of our two main samples, where it is safe 
to assume that the peculiar properties of individual stars 
are likely to average out. However, to determine whether 
both the Li enrichment and enhanced angular momen- 
tum were consistent with the planet accretion hypoth- 
esis, we needed estimates of both the stellar mass and 
radius. The latter can vary by two orders of magnitude 
during RGB evolution. Instead of averaging these prop- 
erties over our stellar sample, we selected three Li-rich 
rapid rotators that are the best examples in our sample of 
stars that may have engulfed planets. Stellar isochrones 
were used to estimate stellar masses and radii. Com- 
bining these stellar properties with the stars' measured 
v sin i and an estimate of the maximum orbital separa- 
tions of accreted planets, we find that minimum planet 
masses of ~ 4.5 Mj up can account for the rotation of the 
two most Li-rich, rapidly rotating stars. This planetary 
mass estimate is comparable to the ~ 6 Mj up needed to 
explain the global Li enrichment seen in the rapid rota- 
tors. The rotation of the third (and least Li-rich) of the 
selected stars can be explained with a minimum planet 
mass of only 1 Mj up . 

However, we also found in the stellar evolution anal- 
ysis that our stellar sample did not reproduce the de- 
tailed -A(Li) and 12 C/ 13 C abundance patterns we ex- 
pected. First, there is significant scatter in the ^4(Li) 
abundances in every evolution group, and we found stars 
with low 12 C/ 13 C at stages earlier than the completion 
of FDU. This scatter may be due to the wide range of 
stellar metallicities and masses represented in our sam- 



ple combined with both a temperature-dependent A(Li) 
sensitivity and a difficulty in measuring precise values for 
12 C/ 13 C > 20. Second, we did not find any classical Li- 
rich stars (-A(Li) > 1.5 dex) at the luminosity bump. The 
most Li-rich slow rotator had A(Li)= 1.12 (and may be 
associated with the luminosity bump, the fourth of our 
six evolutionary classes) . Furthermore, all stars that are 
more Li rich than this are in the second of our defined 
classes — likely at pre-bump evolutionary s tages . This 
finding is reminiscent of the iKumar et al.l (|2011l ) study 
that also found Li-rich giants at pre-bump evolutionary 
stages. Based only on this concentration in T c g, they sug- 
gested that the stars in their sample are red clump stars 
and that Li regeneration may occur during the He flash. 
However, most of the warm (T g< 4600 K), pre-bump Li- 
rich stars in the Kll study had 12 C/ 13 C< 16, whereas 
our Li-rich pre-bump stars generally have 12 C/ 13 C~ 20. 
On the other hand, we noted that there is a slow rota- 
tor in our sample with nearly identical T e g, logg, and 
[Fe/H] as the most Li-rich star in our sample. That slow 
rotator has a low 12 C/ 13 C suggestive of the enhanced 
extra mixing that should only occur at post-luminosity 
bump stages. Together, these two stars appear to be in 
adjacent stages of the Li- regeneration phenomena that is 
thought to occur at they luminosity bump except that 
they are both too hot to be luminosity bump stars. 

In conclusion, the A(Li) and 12 C/ 13 C of our sample 
showed far greater complexity than we anticipated. The 
tendency for the rapid rotators to show Li-enrichment 
implies that either planet accretion or some sort of ro- 
tational mixing has taken place in these stars. If planet 
accretion is not responsible for the Li enriched stars, then 
the fact that these stars are hotter than the luminosity 
bump presents a problem for the Li-regeneration models 
that generally require the removal of the mean molecular 
weight barrier. Many of the slow rotators in our sample 
showed lower than expected A(Li), suggesting that a va- 
riety in either the initial Li abundances or amount of Li 
destruction on the MS may exist. 
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APPENDIX 

The derivation of Equations ^ and (0) is conceptually straightforward. Both A(Li) and 12 C/ 13 C are ratios of the 
number of atoms. In the former, it is the ratio of Li atoms to H atoms. In the latter, it is the ratio of the isotope 12 C 
to the isotope 13 C. 

LITHIUM ENHANCEMENT 

After a planet is accreted, the new A(Li) — defined in Equation (|A1[) — can be calculated by summing the contributions 
of Li and H from the planet and from the stellar envelope, as shown in Equation (|A2[) . 

A(U) new = log(iV(Li) new /iV(H) new ) + 12.00 (Al) 

A(Li) ncw = log((iV(Li)* + N(U) P )/(N(K)+ + N(R) P )) + 12.00, (A2) 

where iV(Li) and -ZV(H) are the number of Li and H atoms, respectively, and the subscripts "new," and "p" refer 
to the star after planet accretion, the star before planet accretion, and the original planet, respectively. Rewriting the 
definition of A (Li) gives an equation for the number of Li atoms as 

AT (Li) = Af(H)10 A(Li) 10- 12 . (A3) 

Therefore, Equation (|A2I) becomes 

A(Li) ncw = log((7V(H)aO A(Li) *10- 12 + ^(H) p 10 A(Li M0- 12 )/(iV(H), + N(R) P )) + 12.00, (A4) 
Rewriting Equation (IA4[) as a sum of logarithms (after factoring out the 10 -12 ) yields 

A(Li) ncw = log(Ar(H)a0 A(Li) * + N(U) p 10 A ^) + log 10~ 12 - log(JV(H)* + N(B) P ) + 12.00, (A5) 
which simplifies to 

A(Li) now = log(iV(H)a0 A(Li) * + N(R) P 10 A ^) - log(i\T(H)* + 2V(H) P ). (A6) 
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The actual number of H atoms comes from the hydrogen mass fraction of the object, A, which is generally known 
(or assumed). That is, Af(H) = XM /ma, where M is the mass of the object, and tor is the mass of a single H atom. 
Assuming that the stellar envelope and planet have the same mass fraction of H (i.e., X* — X p ), then 7V(H) p can be 
written in terms of A(n)* as 

A(H) p = M p /M cnv N(R), = q c N(R)+, (A7) 

where q c = M p /M cnv . Substituting Equation (|A7[) into Equation (|A6|) . factoring out A(H)*, and rewriting as a sum 
of logarithms yields 

A(Li) ncw = log A(H)* + log(10 A(Li) * + q c 10 A(u) -) - logAT(H)* - log(l + q c ), (A8) 
which simplifies to Equation ^ in the paper 

A(Li) now = \og(q c 10 A ^ + 10 A ^*) - log(l + q c ) (A9) 

CARBON RATIO ENHANCEMENT 
In the following equations, we simplify the notation by dropping the N() notation when referring to the isotopes. 
In other words, instead of A( 12 C), we simply write 12 C. After planet accretion, the new 12 C/ 13 C comes from the sum 
of the contributions of each isotope from the stellar envelope and the planet. Thus, 

( 12 C/ 13 C) now = jjgg +1' I (A10) 

For both the stellar envelope and the bulk planet composition, the following relationships are true. We assume that 
all of the carbon is in the two most abundant isotopes. Therefore, 

12 C+ 13 C = N(C). (All) 

If we use r to represent 12 C/ 13 C then Eq. (|A11[) can be expressed in order to solve for 12 C as 

12 C = rA(C)/(l+r) (A12) 

Similarly, Eq. (|A11[) can be expressed as a solution for 13 C as 

13 C = JV(C)/(l + r) (A13) 

The total number of carbon atoms, N(C), in the object comes from the abundance A(C) (the number ratio of C to 
H atoms) and the total number of H atoms from N(H) = XM/m-n: 

N(C) = 10 A(c) - 12 AM/m H (A14) 
Using Equations (|A12|) — (| A14[l with subscripts for the planet and star, Equation (|A10[) expands to become 

r p 10^ c )pl0- 12 ^#^ + raO A(c) *10- 12 X %^\ 

/12p /13p\ _ J_ m H (l+r p ) m H (l+r.) / * T c\ 

( ' jncw " i A(o) pl0 -ia^^ + ^(0)^0-12^^ ■ {Alb) 

Assuming X p w A"^, substituting q c M cnv for M p , and simplifying yields the following expression for the post-planet 
accretion 12 C/ 13 C: 

10 A(C) P iV^ + 10 A ( C ) 



